Commentary 


MANPOWER AND CULTURE 


MANPOWER used to be measured by counting 
reads. Nowadays it is important what is inside the 
heads. Napoleon is supposed to have said, *‘ Only 
aumbers annihilate ’, and to have invented the mass 
ievy though it is doubtful if the concept or the 
practice would have been regarded as novel by 
Genghiz Khan. In the last great war we were told 
that it required eleven men behind the fighting soldier 
to keep him fighting. In times of war there is an 
unanimity of immediate purpose which we rarely 
achieve in peacetime, partly because the immediate 
purpose is paramount. The parameters thus simply 
and starkly laid bare in war, remain coercive in 
peace. It is a dismal way to assess human affairs in 
terms of conflict; we are effete if we pretend these 
factors do not exist. 


A culture is a cultivating. It is more than the 
climate, it is the wHole organic process, spiritual and 
material, aesthetic and functional. In the last two 


hundred years we have developed two cultures which 
appear to some degree to have become polarized, the 
humanist and the scientific, and instead of becoming 
integrated, the two are becoming more polarized. 


Sir Herbert Read, as President of the Society for 
Education through Art is reported to have said’: 
‘The ideal of technology is complete automation; 
the corresponding ideal of education is a human brain 
that controls itself, free from emotional or idealistic 
entanglements, free, above all, from originality of any 
kind. Functional thinking is like functional machinery, 
cold, precise, imageless, repetitive, bloodless, nerveless, 
dead .. . For acentury or two we may be able to control 
its devastating effects by establishing a universal 
tyranny more absolute than any yet experienced by 
mankind. But disaster, as the scientists themselves 
have predicted, is now a mathematical certainty.’ 
Such despair is so wildly ludicrous to anyone 
with any appreciation of the scientific culture that 
it is brushed aside as incomprehensible. 


This despair, this negative preciosity should not 
be brushed aside. It is a widely held view and it is 
thoroughly alarming for the circumstances of which 
Read complains are the very circumstances of his 
opportunity. People who are terrified of survival 
have little time for the enriching of life. The chance, 
the opportunity, is the leisure that power over 
natural resources brings. The population of the 
British Isles was 19 million in 1815 and of the years 
that followed Ramsay Muir writes*: 

‘It is the most poignant commentary on the con- 
dition of Britain that, according to the state of trade, 
from one tenth to one fourth of the population of 
England and Wales were paupers, drawing allowances 
from the poor-rates to supplement their wages because 


these were insufficient to maintain life. What was 
worst of all, by the receipt of poor-relief, men were 
forced to submit themselves and their families to a 
sort of slavery to the poor law authorities . . . British 
liberty had become a very unreal thing for those who 
were subjected to such conditions.’ 
We have many many problems to occupy our 
attention today but there have been some changes. 
Dr MARK ABRaMs has recently prepared a survey 
of teenagers’ spending*; the average British youth 
has £5 a week to spend and his girl £3. They form 
a very important market. 


The purpose of education is to widen horizons 
and an educational system is a reflection and an 
expression of a culture. 

‘The curious thing was that in Germany, in the 
1830s and 1840s, long before serious industrialization 
had started there, it was possible to get a good uni- 
versity education in applied science, better than 
anything England or the United States could offer 
for a couple of generations. I do not begin to under- 
stand this; it does not make social sense: but it was 
so, with the result that Ludwig Mond, the son of a 
court purveyor, went to Heidelberg and learnt some 
sound applied chemistry. Siemens, a Russian signals 
officer, at military academy and university went 
throvgh what for their time were excellent courses 
in electrical engineering. Then they came to England, 
met no competition at all, brought in other educated 
Germans and made fortunes exactly as though they 
were dealing with a rich, illiterate colonial territory*.’ 
Sir Charles Snow has called his Rede Lecture* 

‘The Two Cultures and the Scientific Revolution’. 
Certainly the two have much to learn from each 
other. At present through ignorant jealousy each 
copies the worst features of the other. There is a 
pride in incomprehension which is simply disgrace- 
ful—and stems from laziness and fear. Soon a 
competent journeyman will have the leisure of a 
landed aristocrat and will need the equivalent 
training of a University General Degree. Without 
the technical training he will beilliterate; if his leisure 
is not developed, he will be ‘as the beasts that 
perish’. Our present educational system is a 
patchwork of legacies and pressures. The design of a 
patchwork has the real merit of flexibility; there is 
much that is good but there is a lot of lumber. What 
is needed is a radical re-thinking of the purposes, 
and this is really the problem for there is no single 
purpose. The purposes are as diverse as the culture 
it should reflect. 

' See Daily Telegraph (23rd April, 1957); 

* Muir, R. Short History of the British Commonwealth, p. 314 

London: Phillip, 1934 
® Manchester Guardian (\st August, 1959); 


* Snow, C. P. The Two Cultures and The 


: Revolution 
London: Cambridge University Press, 1959 
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THE USE OF CERAMIC MATERIALS 
FOR GAS TURBINE BLADES 


J. M. PURI 


Imperial Chemical Industries (India) Ltd, Calcutta 


The article is a survey of the use of ceramic materials for gas turbine blades and 

is based partly on the author’s own work at Imperial College, London, and partly 

on a paper presented by T. H. BLAKELEY and R. F. DARLING to the North East 

Coast Institution of Engineers and Shipbuilders. This paper was one of the 
entries in the 1958 Waverley Gold Medal Essay Competition. 


THE outstanding need to raise the operating 
temperatures, and thereby the thermal efficiency, of 
the gas turbine has given tremendous impetus to 
the search for satisfactory materials of construction 
for fabricating blades and other component parts. 
For a number of years work has been done on the 
development of gas turbines employing inlet tem- 
peratures of 1200°C or more. These temperatures are 
well beyond what could be withstood by metallic 
materials, and hence attention has to be turned to 
non-metals, i.e., ceramics. Furthermore the possi- 
bility that certain types of ceramics might meet this 
need has focussed attention on the general question 
of ceramics as engineering materials. 


At temperatures of the order of 1000°C, metals 
often fail as materials of construction through low 
strength, low creep resistance or poor resistance to 
oxidation and, more generally, their pronounced 
instability to other various forms of atmospheric 
attack. In comparison with metals, ceramics are 
immediately attractive, because of their freedom 
from oxidation and their chemical stability at high 
temperature. Their low density (seldom exceeding 
four) is also an attractive feature to minimize centri- 
fugal stresses in parts rotating at high speed. 


Properties and Limitations of Ceramics 


Besides the above mentioned advantages of ceramics 
various other properties may be considered. All 
common commercial refractories and ceramics are 
based on mixtures of some major and several minor 
components. During the burning process, needed 
to produce the ceramic, eutectics are formed which 
melt to give liquid in addition to the solid crystalline 
phases and when the material is cooled, the liquid 
often persists as glass. The strength of this type of 
ceramic, at both normal and high temperatures, is 
due to the inherent strength of both crystals and the 
glass or liquid bond. As the temperature is raised, 


the properties of the liquid phase soon become 
dominant, and as a result the material changes from 
a state of perfect elasticity to one of plasticity. For 
most commercial materials this transition is effective 
below 1000°C, when the material begins to show 
pronounced creep, which is centered in the liquid 
phase, but probably much modified by the crystalline 
components. 

Whilst the strength of many materials of this type 
is adequate at room temperature, it is soon lost on 
heating and, in addition, creep becomes excessive. 
The corollary is that ceramics of suitable strength 
and creep resistance at high temperature will be 
found only in two types of material: 

(i) multi-component materials forming eutectics of 

very high melting point, and 

(ii) single-component high melting point materials 

giving theoretically no liquid below the melting 

point of the component itself. 

Class (i) consists mainly of the less common materi- 
als with melting points in the range 2000° to 3000°C, 
and much more information on the properties of 
such mixtures is required. On the other hand, there 
is now a good deal of knowledge and experience of 
class (ii) materials, which are simpler to obtain and 
to study. 


Before specific materials are described, further 
generalizations about ceramics should be mentioned. 
Their main drawbacks as engineering materials 
include brittleness, that is the virtual absence of 
ductility, and abrasiveness. In general they can- 
not be machined or built up from stock, although 
the shaping or turning of unfired or lightly fired 
material is sometimes possible with small articles. 
Usually the materials must be fabricated by what 
are, in effect, powder methods and then burned at 
high temperatures’. 


The low ductility is a serious disadvantage in 
several respects. Such materials are liable to break 
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zircon—chromic oxide; 


down as the result of local stress concentrations— 
the ability of metals to yield and distribute high 
localized stresses being notoriously absent. Whilst 
the tensile strength of certain sintered oxides 
appears to be high enough to meet estimated 
stresses under turbine conditions, these high values 
are obtained only when great care is taken in testing 
to ensure axial loading and thus to avoid local stress 
concentration’. In other words, failure in practice 
of intricate shapes such as turbine blades might 
often occur even though laboratory tests indicate 
that the strength of the material is adequate. 
Similarly, the low impact strength is an added 
danger, as materials that are liable to fly to pieces, 
rather than to bend, on collision with particles in 
the gas stream, clearly could not be tolerated. 


Materials Considered 


Numerous ceramic materials have adequate hot 
strength but poor thermal shock resistance. The 
fundamental physical properties that favour good 
thermal shock resistance are: low coefficient of 
expansion, low Young’s modulus and high strength, 
and high thermal conductivity. These criteria pro- 
vide an initial sorting test for possible materials; 
they ignore various textural factors which are less 
easy to assess quantitatively. 


The materials considered of interest were 
beryllia, alumina, zircon, silicon carbide, and 
metal—ceramic bodies (sometimes known as ‘ cer- 
mets ’). Of these, beryllia was not considered in this 
investigation® on account of its toxicity and atten- 
dant complications in manufacture. Before going 
on to describe the actual tests performed, a brief 
account will be given of the materials. 


4/umina—Recrystallized alumina was used as a 
blade material. Its thermal shock resistance was 
considerably improved by incorporating relatively 
coarse material; this technique broke up the con- 
tinuous structure. Figure ] (LEFT) shows two blades, 
the bigger of the two blades was coarse alumina 
(bonded with silica), and the other blade recrystal- 
lized alumina. 
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Figure 1. (Lert) Alumina blades; (CENTRE) hollow and solid blade made from 


and (RIGHT) blades made from IC silicon carbide 


Zircon—The properties of zircon (zirconium sili- 
cate) indicate that it is a promising material. 
Preliminary tests on zircon showed that it was 
relatively weak at 700° to 800°C, as compared with 
higher or lower temperatures. It was found that 
this defect could be eliminated by small additions 
of chromic oxide Cr,O, with an optimum at about 
two and a half per cent addition (see Figure /). 

An unfortunate property of zircon is its tendency 
to dissociation into zirconia and silica at, or near, 
its firing temperature. It had been found that 
different batches of zircon bars apparently fired 
under identical conditions had dissociated to differ- 
ent extents. The evidence for correlating creep 
resistance with zircon dissociation is, of course, 
not conclusive, but some support for the general 
view is available. Re-association can be made to 
occur by refiring the material and holding it at a 
lower temperature. Alternatively dissociation can 


be prevented by firing to a lower temperature. On 
the basis of results obtained, a maximum of 1700°C 
for the firing temperature was chosen. Further tests, 
including creep, indicated that zircon-chromic oxide 
fired to 1700°C will probably satisfy the unexacting 
requirements of a single-stage turbine. 





Figure 2. Typical cracks near blade roots 
Figure 1 (CENTRE) shows a hollow and a solid 
blade made from zircon-chromic oxide. 
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Silicon Carbide—This too is a very promising 
material. Commercial silicon carbide bodies are 
usually made by bonding silicon carbide grain with 
clay. This however, adversely affects the properties 
of the pure material. More recently a new form of 
self-bonded silicon carbide, termed ‘ IC’, has been 
developed*, which is more resistant to oxidation and 
distortion than the clay-bonded carbides and it 
combines good thermal conductivity, good high- 
temperature strength, and good creep resistance. 
Blades made from IC silicon carbide are shown on 
the right in Figure /. 


Metal—Ceramic Bodies—It is obviously desirable 
to combine in one material the high temperature 
strength of ceramics with the resistance of metals 
to thermal and mechanical shocks. Such materials 
known as metal—ceramics or ‘cermets’ require 
further investigation. _Cermets such as cobalt- 


bonded titanium carbide have been considered as 
were 


possible blade materials. Tests however, 


Figure 3. Fractured solid blades 
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confined to the metal-metal oxide systems, of which 
a chromium-—debased alumina body was chosen as 
being the one giving the best combination of ease 
of fabrication and desirable physical properties. 


Testing of Blade Materials 


The test rig consisted of a cascade of five blades 
mounted at the outlet of an oil-fired combustion 
chamber which discharged into the throat of a 
water-cooled diffuser. The gas temperature was 
measured by two refractory sheathed platinum 
platinum-rhodium thermocouples just upstream of 
the blades, and there was also provision for measur- 
ing the pressure drop across the cascade. The blade 
cascade could be removed without disturbing the 
rest of the rig. Five blades were cemented in a metal 
channel piece between the refractory holding pieces. 


First Series of Rig Tests 
The first four cascades to be tested in the rig were 
all solid blades, slip-cast in zircon-chromic oxide, 





Figure 4. Fractured hollow blades 



























Figure 5. Blade segment partly 
dismantled 





with the roots gripped firmly between the packing 
pieces. The method of test was to light up as gently 
as possible and slowly increase the air and fuel flows 
until a temperature of 1205°C and a stress of 
1000 Ib/sq. in. were reached. Thereafter conditions 
would be held steady for several hours, or in certain 
cases a thermal shock was administered by abruptly 
reducing the fuel so that the temperature fell to 
about 705°C. On some occasions the stress was 


increased to 2000 Ib/sq. in., maintaining the same 
temperature of 1205 C. All the blades fractured 
after a few hours running under such conditions, 
in several cases without any deliberate thermal shock 
being administered. The fractures were usually close 
to the root; typical instances of such fractures are 
shown in Figures 2 and 3. 
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Figure 6. Blade segment assembled 


It can be shown theoretically that the stresses set 
up in a hollow blade under conditions of changing 
gas temperature are proportional to the wall thick- 
ness, the wall thickness of a solid blade for purposes 
of comparison being taken as half the maximum 
thickness. The blades were fracturing under aero- 
dynamic stresses much lower than the known 
strength of the material, and it seemed that the 
failure was due to thermally induced stresses. The 
next blades to be tested were therefore made hollow 
with a wall thickness of about one sixteenth of an 
inch, but they too failed even more rapidly than the 
solid ones. Figure 4 shows a cascade of such blades, 
three fractured blades and two relatively intact 

On further consideration it became evident that 
the exposed portion of the blades would heat up 
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more rapidly than the root, thus giving rise to 
exceptionally heavy stresses at the junction of the 
two portions, and this effect would be accentuated 
in the case of hollow blades. Therefore when a 
clearance was established, of the order of about 
0-025 in., between the roots and the packing pieces, 
the first failure occurred after six and a half hours. 
After a total of about 19 hours running, two blades 
still survived the tests, which included eight starts 
and stops and four deliberate thermal shocks from 
1205°C to about 705°C. This was a great improve- 
ment, and showed the need for supporting the blades 
loosely between the holding pieces. 


Second Series of Rig Tests 

This series consisted of tests on an experimental 
single-stage turbine with liquid-cooled rotor and 
refractory stator blades, to run at 1205°C. Being 
a single-stage machine, the fixed blades could be 
supported at both ends, thus greatly reducing the 
stresses and generally simplifying the problem of 
blade mounting. The mounting that was devised is 
shown in Figures 5 and 6. The blades are of prismatic 
form and both ends fit loosely into suitably shaped 
recesses in the outer and inner holding pieces. 
Twelve of them fit together to form a complete ring. 


The next step was to test cascades of such blades 
in the rig. The water cooled diffuser could not be 
used in this case, and the test section had to be of 
segmented shape. The blade cascade was arranged 
to exhaust to atmosphere through a simple refrac- 
tory-lined duct, and the rig, with the exhaust duct 
removed, is shown in Figure 7. 
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Modified Test Rig 


The test procedure was much the same as before. 
Owing to the reduced velocity and pressure of the 
gas, the heat transfer coefficients under thermal 
shock conditions were considerably lower than 
would occur in an engine. To make up for this 
temperature difference, the normal practice was to 
reduce abruptly the gas temperature from 1205°C 
to about 705°C. 


The first four or five cascades of blades to be 
tested were all of zircon-chromic oxide. The best 
of these blades showed themselves capable of with- 
standing numerous thermal shocks without visible 
damage. Others, developed cracks, usually in a 
longitudinal direction, as shown in Figure 8. It was 
discovered that the properties of. this material were 
dependent to a considerable degree on the exact 
nature of the firing treatment. Once the correct 
firing conditions were determined and arrangements 
made to control them strictly, blades could be made 
that were apparently unaffected by twelve or more 
thermal shocks. Several such blades were run for 
prolonged periods at a steady temperature of 
1205°C, and after periods varying from 80 to 150 
hours they began to yield to the stress and eventually 
developed cracks as shown in Figure 9. 

All these failures were considered to be due to 
what was essentially a form of creep. This was borne 
out by laboratory tests. 


During the tests described above, the holding 
pieces had suffered much more severely than the 
blades themselves. 


The webs between the recesses 


Figure 7. Test rig 
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Figure 8. Longitudinal cracks 
after 19 hours running and ten 
thermal shocks 





cracked very easily and the overhung portion of the 
inner holding piece tended to break off. These hold- 
ing pieces had been made of zircon, and it was 
decided to change to one of the proprietary mater- 
ials* which was more resistant to thermal shock. At 
the same time a reconsideration of the design 
showed that the pitch/chord ratio of the nozzle 
blades could be increased somewhat without sacrifice 
of efficiency, and the total number of blades was 
therefore reduced from 72 to 60, i.e., five blades per 
holding piece instead of six. The depth of the 
recesses was also reduced, and the combined effect 
of these modifications was to give a more robust 
holding piece with thicker webs. The new type of 
holding blade gave no trouble throughout the re- 
mainder of the rig tests. 


A further material to be tried was a metal—ceramic 
body consisting of 75 per cent chromium and 25 
per cent debased alumina. These blades showed 
themselves entirely resistant to the worst thermal 





Figure 9. Zircon-chromic oxide 
blades showing distortion due to 
creep 


shock that could be applied, but very weak in creep. 
At 1205°C they distorted appreciably in two or three 
hours, and even at 750°C they showed appreciable 
distortion in 15 or 20 hours, as will be seen from 
Figure 10. 


Several self-bonded silicon carbide blades were 
tested and these proved themselves immune to the 
worst thermal shock that could be administered. 
One set was given sixteen shocks without any sign 
of damage appearing, and the test was subsequently 
carried on at a steady temperature of 1205°C for a 
total running time of 363 hours. The only damage 
suffered by these blades was a certain roughening 
of the concave faces, as shown in Figure //. Analy- 
sis of this glassy material, which had evidently been 
molten at the working temperature, showed it to 
contain elements not present in the original blade, 
and it is believed to have come from inorganic 
material present in the combustion gases, not 
necessarily from the fuel. 
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Turbine Tests 


The first of blades to be tested were the zircon- 
chromic oxides ones. It was considered advisable to 
test them first without the rotor in position and the 
turbine was therefore assembled in this condition. 
Figure 12 shows the nozzle blade ring assembled in 
the turbine casing. The complete nozzle ring was 
then run for about 43 hours, including 11 hours at 
maximum temperature. Two of the blades were 
found to be cracked after a short period of running 
and two more after the completion of the above 
period. The remainder of the blades were in excel- 
lent condition. 


Figure 10. Chromium 
alumina blades showing 
distortion 
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The most serious damage was to the holding 
pieces, particularly the inner ones which are not 
supported on the downstream side. The materia! 
originally used was chosen entirely for its thermal! 
shock resistance and it was evident that something 
with a greater strength would be required. 


Conclusions 


It is too early to make dogmatic statements on the 
outcome of the investigation into the use of non- 
metallic materials for stator blades. It is probable, 
however, that the development of metallic materials 
with good strength at high temperature is approach- 





Figure 11. Silicon carbide 
blades after running for 
363 hours 
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Figure 12. Blades assembled in turbine casing 


ing its limit, and if gas turbine operating tempera- 
tures are to be raised to the same thermal efficiency 
level as a diesel engine, alternative types of blading 
must be sought. The choice lies between metallic 
blades, cooled in some manner or other, or non- 
metallic blades operating at the same temperature 
as the gas. 


In turbine rotors the centrifugal field imposes 
high tensile stresses on the blades but also enables 
them to be cooled conveniently by the thermosiphon 
method, and in this case, therefore, the choice is 
likely to be in favour of cooled metal blades. Things 
are far otherwise in the fixed blading, however, where 
the stresses can be kept a good deal lower. Intensive 
investigation, has revealed the difficulty and compli- 
cation of cooling a ring of nozzle blades effectively, 
and this is an obvious field for the application of 
non-metallic materials. There is an inherent diffi- 
culty in the relatively high brittleness of these 
materials, and this must be faced. An adequate 
factor of safety must be built with the design, and 
the designers, long accustomed to metals, must adapt 
their thinking to utilize the available materials, 
which at the best are somewhat akin to cast iron 
(at room temperature) in their behaviour. 


The work described in this essay has indicated 
that stator blades can be made in a ceramic material 
(zircon) but that their factor of safety is low. Metal 
ceramic blades are more promising, but there is a nice 
balance between good creep resistance and good 
high temperature strength on the one hand, and 
good thermal shock resistance on the other. 


More promising are the other materials already 
discussed, such as the newer forms of silicon carbide. 
These have excellent high temperature strength, good 
resistance to oxidation, thermal shock and creep. 
Fabrication of these bodies still presents difficulties ; 
but here also, further improvement can be expected. 


In tests to date the holding pieces have given, if 
anything, more trouble than the blades themselves, 
and the prevention of gas leakage is another serious 
problem which must be solved if progress is to 
continue. The answer here may lie in the use of 
metallic mountings—possibly protected with a 
coating of refractory—which could be relatively 
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easily cooled as only a small amount of surface is 
exposed to the hot gas stream. 


Although much work remains to be done, a promis- 
ing start has been made and as the incentive is great, 
the future can be faced with considerable confidence. 


The author wishes to thank the Editor of the 
Transactions of the North East Coast Institution of 
Engineers and Shipbuilders® for permission to 
reproduce the illustrations with the exception of 
Figure 1. The author acknowledges gratefully the 
assistance rendered to this work by Mr E. S. White, 
of The Morgan Crucible Co. Ltd, and to Mr L. R. 
Barrett, of the Imperial College of Science and 
Technology for his guidance. 
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THE DESIGN, MANUFACTURE AND TESTING 
OF INFRA-RED PHOTOCELLS 


D. H. ROBERTS 


The Plessey Co. Ltd, Caswell Research Laboratories, Towcester 


This is the ninth article in the series on the technology of semiconductors which 
has been published at monthly intervals from October 1958 to May 1959. The 


present article considers photoelectric devices. 


It will be followed by three 


others dealing with other types of devices, intermetallic compounds with 
particular reference to bismuth telluride and organic semiconductors. 


THE INVESTIGATION of the photosensitive properties 
of semiconductors preceded the discovery of the 
transistor by a great many years. Edmonde 
Becquerel discovered the photovoltaic effect in 
1839, and Willoughby Smith was _ investigating 
photoconductivity in selenium in 1873. The first 
thallium sulphide cell was made by T. W. Case in 
1920, and following on from that came the develop- 
ment of lead sulphide infra-red detectors in the 
early stages of the second world war—the most 
rapid development taking place in Germany. 
However, apart from the work of Gudden and 
Pohl in the 1920s and 1930s most of this work was 
of an applied nature, and for a complete under- 
standing of the photovoltaic and photoconductive 
processes it has been necessary to await the stimulus 
given to solid state research by the discovery of the 
transistor at the Bell Telephone Laboratories in 
1948. 


Since that date there has been considerable work 
on the photo-effects in single crystal bodies of 
materials such as silicon, germanium, and indium 
antimonide, and it has been gratifying to find with 
these materials that we can accurately predict the 
observed behaviour—in complete contrast to the 
case of evaporated films of the lead sulphide group. 


There are four types of quantum detectors—as 
opposed to energy detectors such as bolometers— 
but one of these, the photoemissive cathode does not 
strictly rely on the semiconducting nature of the 
cathode, so it will not be considered any further. 
The other three types are photoconductive, photo- 
voltaic and photoelectromagnetic. These three all 
require quanta of energy greater than a certain 
threshold, this being the energy required to produce 
an electron-hole pair, or a single carrier in the 
extrinsic case. The subsequent behaviour of these 
excess carriers depends on the mode of operation. 


Operation of Photocells 


In the class of detectors known as photoconductors, 
the excess carriers produced by the incident radiation 
are observed as an increase in the conductivity of 
the specimen. It can be shown! that the signal 
voltage V, from an intrinsic photoconductor is 
given by* 

y V.Oot (1 e°2) 

* nd (1+ 0%?) 


where V, is the bias voltage, OQ) the mean quantum 
flux density at the surface, t the decay constant of 
excess carriers and a the absorption coefficient of the 
radiation. Z is the thickness, n the equilibrium 
carrier density and w the angular frequency of 
modulation. (The last is equivalent to 2f, where 
f is the modulation frequency.) This equation 
defines the responsivity, usually given in volts per 
watt. It can be seen that given suitable matching, 
carrier mobility is unimportant but that the res- 
ponsivity is oc t/nZ. In the various manufactur- 
ing techniques this fact together with the need to 
maintain a low value of S, the surface recombination 
velocity, is of paramount importance. 


sed 


The Photovoltaic Detector 


Any minority carriers, produced by the incident 
radiation, that reach the junction by diffusion are 
extracted by the local electrostatic field to become 
majority carriers. In the usual short circuit or 
matched load conditions a current flows through 
the load to preserve charge neutrality in the p- and 
n-regions. (The photovoltaic effect is shown in 
Figure la and /b.) 
The short circuit current i, is given by 

‘ AeQo aula .(2) 

where 7 is the collection efficiency and is a function 


* A list of all the symbols used wil! be found in Appendix J on page 264 
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c’ the diffusion length, the decay constant of the 
excess carriers, the surface recombination velocity 
S, the junction depth and the absorption coefficient 
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The photovoltaic effect 


for the radiation. In regions of high absorption co- 
efficient, 7 is of the order of 0-65 if S is assumed 
negligible and the junction depth equals the diffusion 
length of electrons in p-type material. These are 
then the basic design parameters: to maintain a low 
surface recombination velocity and control the 
junction depth. 











Figure 1. 


Photoelectromagnetic Effect 


The photoelectromagnetic effect (see Figure 2) is the 
most recent mode to have been used in a practical 
photocell*. Non-penetrating radiation produces a 
concentration gradient of free carriers and hence a 
diffusion current of electrons and holes from front to 
back. If a magnetic field is now applied perpendicular 
to the plane of this diffusion, the positive and nega- 
tive charge carriers are deflected in opposite ways, 
thus producing a voltage difference between the ends 
of the semiconductor. 

If one assumes that the surface recombination 
velocity S at front surface is negligible, S at the back 
tends to infinity and the thickness of the specimen 
will be about the same as the diffusion length. One 
then obtains the maximum photocurrent i,,, 

: AQ eB 
Ise = 
where Lp is the electron or hole diffusion length, 


(Lpette — Lpnltn) eee 
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u the electron or hole mobility, A the receiving area 
of the detector, B the magnetic flux and YX is the 
inter-electrode distance. And again, these con- 
ditions set our design parameters, it being necessary 
to control the surface recombination velocity both 
at the back and the front, and specimen thickness 
(apart from control of the material properties, 
lifetime, impurity concentration, carrier mobilities). 


It is worth noting one important difference bet- 
ween the photoelectromagnetic and photoconduc- 
tive modes, namely the dependence on thickness. In 
the first case specimen thickness should be about 
equal to the diffusion depth whereas in the second 
case the thickness should be nearly equal to I/a, 
where a is the absorption coefficient for the 
radiation. (Some of the advantages of the photocon- 
ductive and photoelectromagnetic devices are com- 
pared in Table /.) 


Noise 


Apart from the response of a photocell to a given 
amount of radiation its practical performance 
depends on the noise level in its operating condition. 
There are four types of noise occurring in photocells; 
more complete discussions of noise in semiconduc- 
torsare given by R. E. BurGeEss* and B. L. H. WILSON’*. 


The first, known as the Johnson noise, is produced 
by the voltage fluctuations across any two terminal 
impedance Z in thermal equilibrium at absolute 
temperature 7. At normal frequencies and tempera- 
tures the noise is given by 

Vin = 4kT.JAZ).Af oer 
where k is Boltzmann’s constant, and Af is the 
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band width, and the Johnson noise can obviously 
be reduced by lowering the temperature, the detector 
impedance, and the bandwidth. 


There is excess noise, called Shot noise, in a 
semiconductor carrying a current compared with a 
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Table 1. Comparison of advantages of photoconductors and photoelectromagnetic detectors using indium antimonide 





Photoconductors 


Photoelectromagnetic detectors 





(i) Superior spectral response—edge occurs at 
slightly longer wavelength 


(ii) Less bulk and weight since no magnet is re- 
quired 





(i) Requires no bias supply—important in tran- 
sistorized equipment 


(ii) Easier to manufacture because of the thickness 
effect discussed earlier 


Likely to yield the more sensitive detector— 
because of greater responsivity of the photo- 
electromagnetic effect and lower noise level due_ | 
to absence of bias current 


(iv) Less dependent on lifetime and surface recom- 
bination velocity 


(v) Easier to prepare suitable starting material, | 
p-type doped rather than high purity intrinsic | 
being used 


(iii) 





metal and this is due to the fact that the excitation 
and recombination processes governing the number 
of free charge carriers are randomly occurring dis- 
crete events. (Burgess*® gives a full treatment of 
shot noise.) Provided the lifetime of excess carriers 
is very much less than the transit time, it can be 
shown that for an intrinsic semiconductor, shot 
noise can be described by 

ey 2el.(t/t'). Af eye 
where / is the d.c. current bias, t the decay constant 
of excess carriers, t’, the transit time and Af the 
bandwidth. 


It is found that there exists a further type of 
current dependent noise, known as flicker noise. 
The origin and detailed mechanism of this type of 
noise are very much open to speculation®. Flicker 
noise may be described by the equation 


E'l"f-“ Af er), 


j2 
TEN 


where uw and w are constants and approximately ~ 


equal to two and one respectively. Some data on 
the dependence of E’, and u on the properties of the 
surface and the presence of intercrystalline potential 
barriers was given in a paper by D. H. Roperts and 
B. L. H. WILSON at the symposium on Noise in 
Fixed Resistors, London, March 1958. 


Photon noise is in a way analogous to shot noise, 
since it is due to the quantum nature of the incident 
radiation. In a photoconductor 

2, = 2Cel Af sna 
where /, is the d.c. photocurrent produced by Qo 
and C, is a constant. Photon noise in the signal 
is the ultimate noise in any detector, but, in 
room temperature detectors this is unlikely to 
be observed compared with the other three types. 
Photon noise is also generated by the ambient 
radiation particularly important in long wave- 
length detectors. The increase in sensitivity of 
lead sulphide and lead telluride cells on cooling 


their surroundings has been given by T. S. Moss’ 
as evidence that they are photon noise limited. 


Manufacturing Techniques: 
Photoconductive 


Very sensitive visible radiation detectors made of 
cadmium sulphide have been made using powder 
techniques, with either a painting system with 
polycrystalline powder® or sintering a film of 
amorphous powder’. Infra-red photoconductors 
however have been made principally by sublima- 
tion’®, chemical deposition’! and the use of lapped 
and etched single crystal filaments'*. These will 
now be considered in more detail below. 


Evaporated Film Method 

The evaporated film method has been widely used 
with the lead sulphide group of photoconductors. 
As an example lead selenide films have been pre- 
pared by evaporating starting material of the form 
Pb;Se,03. This material was prepared by mixing 
the appropriate amounts of lead, lead oxide PbO, 
and selenium, sealed in a pyrex tube at a pressure of 
less than 10-5 mm, fired at 450°C and the final 
product was crushed and sieved. A few milligrams 
of this starting material were then placed in a 
pyrex glass cell blank, with bubble window and 
aquadag-tungsten electrodes, which was sealed to a 
vacuum system and evacuated to better than 10-> mm 
of mercury. The actual sublimation process was 
performed by flame-heating the base of the cell 
blank, causing a film of lead selenide to be deposited 
on the window between the electrodes. 


This lead selenide film then required sensitizing 
by heat treatment in air, the choice of time, tem- 
perature, and air pressure being determined by the 
characteristics of the film as prepared and the 
whims of the operator! During the sensitization, 


the conductance and photoconductivity of the film 
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uader standard illumination would be measured in 
a simple d.c. circuit. The variation of these and 
other parameters, thermoelectric power and carrier 
lifetime, with sensitizing treatment is given in a 
recent paper’. At the appropriate stage the cell 
blank would be evacuated and sealed off ready 
for use. 


: Indium arsenide filament 
Copper wires 


™~ . 
Insulating 
Copper layer 


block 











Figure 3. Indium arsenide photoconductive cell 


Chemical Deposition 

Chemical deposition methods have been used and 
again lead selenide will be used to give a specific 
example, although this method has been most 
widely used on lead sulphide by B.T.H. Ltd in this 
country and Eastman Kodak, among others, in 
the United States. Films of lead selenide are repro- 
duced from the colloidal phase in the reaction be- 
tween aqueous solutions of a lead salt and seleno- 
urea in the presence of a base. A glass substrate 
having been degreased and washed was rotated in a 
solution containing lead acetate and seleno-urea, 
hydrazine hydrate being used as the base additive 
to control the reaction!!. Having deposited a film 


of the required thickness the substrate was removed - 


from the solution, washed, and dried. A suitable 
electrode system could then be painted on the film 
using aquadag or acolloidal gold suspension. Workers 
using lead sulphide have subjected their chemically 
deposited layers to a tempering or sensitizing treat- 
ment by baking the deposited films in air at atmo- 
spheric pressure—sometimes for quite long periods. 


Solid Filaments Technique 

Whatever failings the solid filaments technique 
may have, it does possess one major virtue. It has 
been possible with materials such as indium anti- 
monide!? and indium arsenide to construct photocells 
whose performance could be understood and whose 
sensitivity could be predicted theoretically with a 
fair degree of accuracy. The method described below 
is the most recent application of this technique to 
indium arsenide. 


A zone refined ingot of indium arsenide was cut 
into slices, these slices being lapped to a thickness 
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of about 0-1 mm and then diced to give filaments 
3x1 mm which were suitable for mounting and 
etching. In order that the detector shall not be 
dissipation limited! it is necessary to provide an 
efficient heat sink. This was done as shown in 
Figure 3 by mounting the indium arsenide filament 
on a copper block with a thin insulating layer 
isolating the filament electrically whilst permitting 
good thermal contact. Wire contacts were soldered 
on to the filament to connect it across a co-axial 
socket. After mounting, the region around the 
contacts was masked with a hard setting resin or 
wax and then the specimen was etched to reduce 
the thickness toward the ideal of a few microns and 
to reduce the surface recombination velocity. 


Measurements performed on this type of detector 
indicated that in the near future it will compete 
with lead sulphide cells for use in the 1 to 3:5 
micron region. 


Manufacturing Techniques: 
Photovoltaic Detector 


The photovoltaic detector has certain virtues (and 
disadvantages) compared with the photoconductor 
which are briefly summarized below. 
(i) The photovoltaic effect is less dependent on 
lifetime and surface recombination velocity. 
(ii) It is not necessary to reduce the overall thick- 
ness to such mechanically difficult levels. 
(iii) The devices are not so easily de-sensitized by 
an increase in temperature or background illum- 
ination. 


Incident radiation 
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Figure 4. Silcon photovoltaic detector 


The disadvantages are two fold: more advanced 
materials technology is required to produce a p-n 
junction and for long wavelength detectors, ¢.z., 
indium antimonide, it is necessary to dope so heavily 
to obtain a p—n junction at room temperature that 
the carrier lifetime and mobility are degraded very 
considerably. 


As an example, the method of making a silicon 
photocell will be described. This is basically the 
same technique as is used to construct a solar battery 
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for the conversion of the sun’s radiation into 
electrical energy. Assuming one started with a 
single crystal of p-type silicon it is necessary to 
diffuse an n-type impurity, i.e., phosphorus into a 





Figure 5. Infra-red detector (by courtesy 
of Technical Ceramics Limited) 


slice of silicon to obtain a p—n junction. In order to 
make a device however it is necessary to make 
electrical contacts to both the p- and n-regions. A 
very convenient way of doing this is to use what is 
known as the electroless nickel plating technique, 
but in practice it is difficult to make an ohmic con- 
tact to high resistivity silicon, so that a p**+*+ region 
is produced in the base region on to which the 
nickel can be plated. 


The procedure, therefore, is as follows. Boron is 
diffused into the slice'* and the boron doped region 
is removed from one side by grinding. Phosphorous 
is then diffused into the slice and the resulting oxide 
film on the surface removed by hydrofluoric acid. 
After nickel plating the whole slice the active regions 
were exposed by etching in nitric acid, and the 
junction was etched in a mixture containing equal 
parts of the two acids, having masked the remainder 
of the structure. After washing and degreasing the 
nickel plating was tinned, using a dip soldering 
technique; nickel tape leads were then soldered on 
and the cell potted in a hard-setting resin. The final 
structure is illustrated in Figure 4. 


Manufacturing Techniques: 
Photoelectromagnetic Detector 


The methods used for the manufacture of photo- 
electromagnetic detectors have somewhat limited 
application, it being unlikely that it will be used 
with any material other than indium antimonide. A 
comparison between this and the photoconductive 
type of detector is given in Table /. 

As already mentioned the first stage in the 
construction of indium antimonide photoelectro- 
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magnetic detector is to cut and lap the single crystal 
to produce filaments 3 x 1 x 0-1 mm. ‘These are 
mounted on an insulator with contacts attached, as 
there is in this case no need to consider dissipation. 
Once mounted the specimen was etched in, for 
example, a mixture of hydrofluoric acid and hydro- 
gen peroxide. In early development the procedure 
was successively to etch and measure the perform- 
ance until the optimum thickness is reached and 
once a crystal had been characterized a large number 
of detectors could be made by etching for a standard 
time. The specimen was then mounted in a magnet, 
a final assembly being shown in Figure 5. The 
magnets used are made of Alcomax III with mild 
steel yoke and permandur pole pieces shaped to 
give a field of about 10,000 gauss in a gap of 1-5 mm. 


Methods of Testing Photocells 

The spectral response is the variation of responsivity 
with wavelength, and in principle it is only neces- 
sary to measure the ratio of input and output signals 
at various wavelength, the input, i.e., radiation, 
signal usually being measured with a_ thermal 
detector of known responsivity such as a Schwartz 
radiation thermocouple. In practice there are many 
variations and many precautions must be taken. 
For the single beam system one must chop at the 
highest possible frequency to differentiate between 
photoconductive and bolometric effects and also 
chop at the input slit of a monochromator to avoid 
spurious signals. For double beam systems the exact 
method used depends on the equipment available. 
Some of the variables are beam splitting, measure- 
ment of signal ratio, radiation source and area of 
thermal detector. Beam splitting involves combining 
with chopping at the output slit or placing a mirror 
in half the optical field or a half-silvered mirror in 
whole of the optical field. The last of these suffers 
the least disadvantages and is the most reliable. 


Measurement of signal ratio is carried out by 
measuring two signals individually and computing 
the ratio at fixed points, by applying both signals 
(rectified) to a potentiometer recorder connected as 
a ratiometer or by using the reference signal to con- 
trol the source output, or slit width, te give constant 
energy at all wavelengths on the detector under test. 
If the last method is used, a tungsten filament lamp 
with a suitable window to transmit infra-red, can 
be used as a radiation source. Otherwise a Nernst 


filament is suitable provided care is taken to shield 
the filament from draughts. To determine the area 
of the thermal detector care must be taken, in a 
double beam system, to ensure that all the radiation 
signal is incident on the reference detector. This is 
easier if the detector area is comparable with the 
photocell under test. 
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Finally it is worth noting that such a technique 
will give a curve of responsivity in volts per watt 
ersus wavelength, whereas from the _ semi- 
-onductor point of view volts per number of quanta 
‘3 More meaningful and the simple conversion is 
worthwhile. 


Responsivity and Noise 


fo determine the responsivity and noise, and hence 
the equivalent noise power, the signal source used 
is ‘a 200°C black body’. This consists of a tem- 
perature controlled furnace whose interior can be 
seen’ by the detector through a series of baffles 
whose geometrical arrangement is such that the 
final aperture behaves as a perfect black body at 
200°C, thus allowing the radiant flux at a standard 
distance to be calculated quite easily. The radiation 
output is generally modulated at 800 c/s and the sig- 
nal from the detector placed at a standard distance, 
20 cm, is fed into a tuned amplifier with a band- 
width of 50 c/s. The noise and signal levels are then 
measured, in the photoconductive case as a function 
of bias current. The ratio of input and output 
signals, at the detector, gives the responsivity in 
volts per watt. This together with the noise level 
allows the equivalent noise power, defined as mini- 
mum incident energy to give a signal-to-noise ratio 
of unity when measured in unit bandwidth, to be 
calculated. 


With photoconductive cells where flicker noise is 
suspected it is essential to check the current and 
frequency dependence of the noise voltage as these 
will determine the choice of operating conditions. 


Response Time 
It is necessary to have a pulse of radiation or a 
modulated sine wave of variable and high frequency 
to determine the response time. The simplest 
method is probably to use a triggered spark dis- 
charge as source and display the signal on a cathode 
ray tube with calibrated time base, the response 
time, t, being defined as the time for the signal to 
fall to 1/e of its initial value. Care has to be taken 
to guard against spurious time constants. If a fast 
sine wave modulation is used there are two possible 
approaches: to measure the frequency dependence 
of responsivity R 
Ro 
a= (1 + w?r?)t 
or to measure the phase difference 5, between input 
and output signals 
6 = tan-twr eT 

and hence calculate t from (8) or (9) where @ is 
the angular frequency. 

It is worth mentioning that these same methods 
can be used to measure the carrier lifetime in a bulk 


. (8) 


semiconductor, for example, silicon, but it is neces- 
sary to take precautions. In the case of traps the 
effect of d.c. illumination on the shape of the decay 
curve must be observed. The short wavelengths are 
filtered out and again the effect on the decay curve 
due to surface recombination is observed in order 
to measure carrier lifetime. Lastly to measure the 
illumination level, it is necessary to standardize on 
signal level as t can be intensity dependent. 
Conclusions 

In the last few years the range of available infra-red 
detectors has expanded quite considerably, as also 
has the understanding of their performance and its 
dependence on basic material properties. It is to be 
expected that the next few years will see an even 
greater expansion in the application of these photo- 
cells in industry. In particular the use of these 
devices for non-contact temperature measurement 
enables surface temperatures to be measured down 
to only a few degrees above ambient with, for 
example, the indium antimonide detector. Further- 
more the combination of this type of detector with 
multilayer interference filters for the infra-red offers 
the possibility of making cheap, small, robust, 
sensitive instruments for gas analysis and chemical 
process control. 

In a different field, the combination of photocon- 
ductive and electroluminiscent elements is already 
being heralded as opening new paths for the com- 
putor designer and time will possibly bring an 
* optical ’ computor. 

The author is grateful to the Plessey Company 
Ltd for permission to publish this article. 
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ELECTRICITY GENERATION 
FROM MOVING GAS STREAMS 


R. G. VOYSEY 


Chief Scientists Division, Ministry of Power, Thames House South, Millbank, London, S.W.1 


Considerable interest has been shown in thermal-electric generators with no 

moving parts other than a gas stream because it would be possible to operate 

such devices at very high temperatures which could well correspond to high 

thermodynamic efficiencies. The author points out that although the magneto- 

dynamic form of ionic generator shows little promise, an electrodynamic 

generator rather like a van der Graaff machine but using ionized particles has 
potential applications. 


FARADAY’S discovery of the magnetodynamic 
generation of electricity was prompted by Arago’s 
observation of the deflection of a magnetic needle 
near a rotating copper disc. Arago claimed that 
the effect held with a rotating conductor of any sort, 
including liquids and gases. This was a bold claim 
since the resistivity, even of a Bunsen flame, is 
nearly 10" times that of copper. Faraday reports! 
that Babbage and Herschel could not demonstrate 
Arago’s effect with any non-metal except gas-coke 
and he himself could not obtain significant results 
with sulphuric acid or salt solutions. However he 
later obtained erratic readings of the flow of the 
960 ft width of the Thames at Waterloo Bridge by 
measuring the voltage induced by the motion of the 
brackish water through the earth’s field!. 


Since Faraday’s time, and particularly after the 
clear formulation of the conception of ions and 
electrons, repeated patent proposals have been made 
for ionic generators using not metals but fluid, 
preferably gaseous, conductors to secure some 
special advantage of efficiency or applicability. 
For example, 1917 was a year of great interest 
in trans-Atlantic radio transmission and A. YOUNG? 
suggested using a fluid conductor device for the 
generation of low frequency radio waves. In the 
next few years patent applications for ionic genera- 
tors appeared in several European countries?~*. 
These and subsequent inventors had little conception 
of the problems but were attracted by the apparent 
possibility of extremely compact thermal-electric 
generators having no moving parts other than a gas 
stream and susceptible to operation at very high 
temperatures, which might imply a correspondingly 
high thermodynamic efficiency. 

The magnetodynamic form of ionic generator is 
probably unattractive but an electrodynamic form, 
rather like a van de Graaff machine but using 
ionized particles might be feasible if properly 
designed. While cold, van de Graaff type machines 
using dust to carry charge are well developed 


(early investigators were Armstrong and Faraday’) 
but there is little evidence of thermal-electric 
machines with ionic current carriers being made. If 
experimental work has been done it probably exposed 
two serious difficulties. These are the high thermody- 
namiccost of providing sufficient current carriers in the 
gas stream and the possible instability of these ionic 
generators to internal disturbances of various kinds. 
These difficulties are worth examination after con- 
sidering the principles and the simplest forms which 
such devices might take because there is renewed 
interest in efficient light weight thermal-electric 
generators of all sorts and a growing appreciation of 
the subtler points of magneto-hydrodynamics which 
might encourage a practical ionic generator. 


As a device employing no machinery, other than 
its auxiliary compressor or boiler, it invites com- 
parison with the fuel cell, the thermocouple genera- 
tor and the thermo-electron generator. The first 
two have been the subject of articles in RESEARCH*”® 
while the third is being actively studied in the 
United States’®-'*. Unlike these three, the voltage 
output of the ionic generator is not controlled 
internally by considerations of the material work 
functions and, in the electrodynamic form, prefers 
very high voltages controlled by external fields. In 
theory the ionic generator suffers from being a semi- 
mechanical heat engine and not truly a direct 
generator of electricity from heat. Its promise of 
high temperature operation makes this considera- 
tion far less important than practical difficulties. 
If feasible it should be far more compact than the 
other three. 


Principles of Generation 


The particular class of electrodynamic machines 
here called ionic generators collect the charge of 
gaseous ions, i.e., charged carriers, propelled against 
electrostatic or magnetic fields, or combinations of 
the two, by a moving gas. The electrical output is 
maintained by a back pressure on the gas flow. 
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Such devices have been suggested, operating in 
reverse, as ion rockets!® and also, working in reverse 
with liquid metal conductors, they are electro- 
magnetic pumps. 

Before considering physical forms of the device it 
is possible to make an estimate of the required degree 


of ionization of the gas for various voltages of 
veneration. 
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Figure 1. Heat engine cycle—ideal con- 


ditions only and not necessarily achievable 
in practice 

Taking some specific thermodynamic cycle, for 
example the very speculative Joule cycle of Figure /, 
a figure can be chosen for the desired isentropic heat 
drop in the device, in this case 240 C.H.U.s per 
pound. Irrespective of the machine’s size or mass 
flow rate this leads immediately to a relation between 
the necessary ion population and the voltage of 
operation. Table /, based simply on the conserva- 
tion of energy, shows the tolerable scarcity ratio of 
ions, i.e., the number of neutral molecules per ion 
which must not be exceeded to carry an adequate 
current at various voltages for 100 per cent efficiency 
of generation. Even a simply ionized molecule, in- 
volving the detachment of one electron, requires the 
energy equivalent of about 20,000°C per molecule. 
Table 1 shows that on the most optimistic expectation 
of operation, a voltage near 30 kilovolts implies a loss 
of energy in ions at the exhaust equivalent to 20°C 
when averaged over all the molecules. At 600 volts the 
loss would average 1000°C and compare with the 
isentropic temperature drop assumed in the cycle of 
Figure 1. A further loss is entailed in providing 
electrons at the output electrode to neutralize the ar- 
riving positive ions. 

In small scale application efficiency is often less 
important than extremely light weight and simplicity, 
but quite high voltages would still be demanded. 
This requirement almost certainly disposes of the 
simplest form of ionic generator suggested (see for 
example Figure 2) but is less damaging to the device 
suggested in Figure 3. These forms are described 
under the next heading. 


Cc 
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Elementary Forms of lonic Generators 


The conception of Figure 2 is the favourite of 
inventors. The gases in passing through a suitable 
expanding nozzle also travel through a transverse 
magnetic field. Under the combined action of the 
magnetic field and the gas velocity any charged 
particles in the gas will try to traverse orbits having a 
circular motion plus a steady drift in a direction at 
right angles to both the gas flow and the magnetic 
field. They will also suffer repeated collisions. With 
a possible and proper choice of design conditions the 
sideways drift preponderates over the circular 
motion and some carriers are imagined to be 
collected finally by side electrodes. The theoretical 
maximum voltage generated is simply the number of 
magnetic lines cut per second, divided by 10°. Even 
with fields of 80,000 lines per square inch and gas 
velocities of 1000 feet per second, the potential 
generated can only be of the order of ten volts per 
inch width of nozzle. 


A number of considerations limit the nozzle width 
so, if the device is to be efficient, it must produce 
extremely large currents to harness the very high 
rate of energy conversion from heat to velocity which 
takes place in the nozzle. The large current in turn 
implies a high proportion of ions in the gas stream 
so that this idea might only be made to work at all 











Table 1 
Neutral Required Average temperature 
molecules voltage loss with 20,000°C 
per ion ions 
(in kV) (in °C) 
50,000 1515 0-4 
20,000 606 1 
10,000 303 2 
5000 151-5 4 
1000 30-3 20 
Ee 








with elaborate and energy-consuming ion injectors 
and collectors. Without them it would give only a 
tiny fraction of its theoretical output and corres- 
pondingly little resistance to gas passage through the 
nozzle. 


Inventors are tempted to believe that some special 
circumstance, for example the use of some easily 
ionizable medium like caesium vapour can help. 
It is true that at the gas temperatures in which we 
are interested a small amount of alkali metals (or 
better, the alkaline earths) give a disproportionate 
increase in conductivity. This still leaves the con- 
ductivity very low indeed and the gas also contains 
charges of both signs which is useless for this device. 
These matters are discussed further under the head- 
ing concerned with the production of ions. It is 
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worth repeating that the table is based on 
simple arguments of energy balance. The results are 
expressed more generally by the following relation: 
the required fractional concentration of ions is 
approximately 0-00043 WM/V. (W is the specific 
output of the gas in C.H.U.s per pound flow, M 
is the molecular weight of the substance and V is 
the voltage of generation.) 

The second simple conception of an ionic genera- 
tor is the electrodynamic one of Figure 3. Advantage 
is taken in this of the high dielectric strength of 
gases, even at high temperature, to work at very 
high voltage. The ions are propelled against the 
repulsion of the collecting electrode. Modern work 
(see reference 17 and its sources) has shown the 
general design requirements (particularly for the 
geometry and terminal arrangements) of normal van 
de Graaff machines and these need not be treated 
here. The gases in the ionic generator would flow 
at over ten times the highest speed of the charging 
belts in the normal van de Graaff. Elementary con- 
sideration shows that while the gas velocity may be 
high and the distance between starting and collecting 
electrodes can be made so small that the transit time 
might be of the order of ten microseconds, a useful 
charge carrier might—and would desirably—suffer at 
least a thousand accelerating collisions in transit. 
There is then the difficulty that an avalanche of ioni- 
zation might set in with the production of a conduc- 
tive path in the stream through which the electrical 
energy-could arc back. This is part of the general 
problem of instability but before it is discussed 
further it will be advisable to clear away the simpler 
matter of suitable thermodynamic cycles. 


Suitable Thermodynamic Cycles 


The arguments behind Table / do not prejudice the 
choice of medium. If any sort of ionic generator is 
ever made, some practical requirement such as the 
ionizing mechanism might restrict the medium but 
there is no criterion for limitation at the present 
stage. It is therefore possible to consider cycles using 
any condensible or non-condensible gas. As with 
most thermal-electric devices it is helpful in con- 
sidering efficiencies to separate the ideal thermal 
cycle from the electrical process. With an ionic 
generator one can think of the appropriate cycle, 
say Rankine’s for condensible ‘media—water, 
mercury, sodium, caesium efc—or the Joule cycle 
(as in Figure 1) for gases like air, combustion gases, 
monatomic gases efc. With these two cycles the 
generator component is the expander-plus-dynamo 
of a steam or gas turbine cycle and these are well 
studied fields. 


It would be out of place to review the extensive 
literature on these cycles and their available media. 
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Briefly, if practically achievable, the ionic generator 
will derive great benefit from its lack of machinery 
and enjoyment of high temperature. The Joule cycle 
conditions of Figure / suggest a target cycle efficiency 
of 83 per cent. A Rankine cycle using some suit 

able low pressure vapour, say mercury between 750°C 
and 120°C, cansurpass 50 percent efficiency. Itis also 
possible to think of the generator as incorporated in 
the relevant stage of any reciprocating cycle, e.g., in 
the Otto or Diesel cycle, with a further expansion of 
its imaginable embodiments. 


Most inventors or promotors of an idea make the 
mistake of claiming overwhelming advantages which 
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Figure 2. 


Magnetodynamic generator 


may very well be real in the long run but compromise 
the early limited success which might be achieved by 
avoidance of immediate problems of deficient know- 
ledge and technique. The basic thermal cycles which 
an ionic generator might use are so efficient that it is 
advisable to exploit them to relax design difficulties 
elsewhere. Thus the device might be developed as an 
inefficient light weight device or only as a super- 
position unit, drawing a valuable few per cent of 
energy from a hot gas and exhausting to a more 
normal heat engine component like a boiler or 
turbine. 


It is difficult to imagine the device working con- 
sistently on other than ash-free fuels. 


Production of lons and 

the Possible Instability of the Device 

It may have been noticed that the word ‘ plasma’ 
has been avoided so far. It suggests a more extreme 
condition of a conducting fluid in which most of the 
molecules are ionized at least singly, carriers of both 
sign are present, and the gas is in something near 
thermodynamic equilibrium at 15,000° to 20,000°C 
or above. None of these conditions is either necessary 
or desirable here unless the generator is imagined to 
derive its energy from low cost nuclear sources. 
Every re-inventor of the device thinks of radio- 
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activity as an easy source of ionization but it is an 
ill matched source. Ionization needs collision and 
radioactive bombardment represents an extreme 
case of using particles of very many electron volts 
to liberate electrons from bonds of only a few elec- 
tron volts energy. 


At the other extreme inventors have hoped for 
sufficient ionization by gentler collisions at thermal 
velocities. A proportion of ions may be produced 
even in cool gases by chance collisions of excep- 
tionally fast molecules but the concentration could 
not rise to useful values, even for kilovolt operation, 
until the gas temperature was above 2000°C. 


While the curious behaviour of hot gases in a 
magnetic field has been noted for over a century, 
see, for example, Faraday’s experiments’, the 
physical conceptions for explaining their electrical 
behaviour were lacking till this century. Even now, 
while technology spreads towards the application of 
very high temperatures, it is still impossible to 
estimate the condjtions of ionization in hot gases 
except under very restrictive assumptions. The 
first attempt, by Saha for estimating the conditions 
within a star, led to an equation which is discussed 
in a chapter on the ionization of flames’®. 


The Saha equation at the outset assumes equal 
numbers of positive and negative ions and involves 


the product of these numbers to evaluate an Ga 


equilibrium constant for their partial pressures. if 
a vanishing proportion of negative ions is assumed, 
the equation would justify the stability of a very 
large proportion of positive ions. In the ultimate, 
such a cloud would be self-repelling, ignoring any 
‘pinch’ due to motion, but the small electric 
pressure could easily be balanced by potential on 
the walls. A more obvious difficulty might be the 
encouragement of electron emission from both walls. 
It seems probable that this will always take place 
to some degree even if the wall is built up from a 
series of guard rings of graded potential and it 
implies a limitation of working temperature and the 
use of materials of very high work function for every- 
thing except the electrodes. It also discourages the 
scattering of alkalis or alkaline earths into the gas 
stream. It has already been said that the advan- 
tages of such additions are probably not real at the 
required conductivities. The device could not rely 
on thermal ionization alone and indeed the latter 
would, if significant, destroy its stability. The 
feasibility of operation rests essentially on economic- 
ally producing a good supply of ions of one sign 
and then taking advantage of the ‘freezing’ of 
electronic reactions at the comparatively low 


temperatures to complete the few microseconds of 
travel through the nozzle, 
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Considering the simple electrodynamic generator 
the charged particles can only proceed against the 
voltage gradient by replenishing their kinetic energy 
from uncharged neighbours in their decelerating 
path. A number of collisions approaching 1000 
would meet the case without running more than a 
negligible risk of meeting sufficient very high energy 
particles to initiate avalanching. There is sufficient 
design freedom in the problem to keep down to this 
order of number of collisions. 


In view of the restricted prospects of the magnetic 
conception of Figure 2 there is no point in discussing 
the very complicated ion behaviour which should 
take place in it. 


Even if the electrostatic form of Figure 3 is thought 
promising, design is complicated by the interaction 
of the generated current with itself and the surround- 
ings (as with armature reaction problems in rotating 
machinery or ‘ pinch’ efc in plasmas proper). While 
Maxwell’s equations tidy up electromagnetic be- 
haviour very nicely, it is difficult to reconstruct 


sani grid 






Expanding 


+ Exhaust 
channel % _ 


t 
2 





s 
Ne t | 
| 


le——QOutput voltage——+ 


Figure 3. Electrodynamic generator 


physical effects from these laws. We can now follow 
nuclear pioneers and take over their developed ideas 
of ‘ pinch ’, ‘ magnetic mirrors ’, ‘ magnetic bottles ’ 
etc but these effects are very much modified in the 
generator’s region of relatively low temperature, low 
ion density and mobility. Consideration suggests 
that with sensible proportioning of the nozzle there 
should be no fundamentally unsolvable problems of 
loss of ions or pinch. The more obscure electro- 
magnetic effects may be exploitable. The 1958 
Geneva Conference on Nuclear Energy contained 
many papers on the dynamics of plasmas proper 
but only a few of the papers are at all relevant here; 
a good general reference is the symposium*’. 


Conclusions 


Despite repeated attempts to patent an _ ionic 
generator no one appears to have made any practical 
form of the device. This is not surprising in view of 
the daunting complexity of design and operation 
which underlies the superficial simplicity of its 
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conception. Yet its basic thermal efficiency and its 
high power—weight ratio are so promising that it is 
surprising that so little experimental work has been 
done. Some exploratory measurement on ion densit- 
ies was made by M. W. TurING and others at the 
University of Sheffield in 1958 with the support of 
the National Research and Development Corpora- 
tion. 


Much more work could be done on the electro- 
static form which may have early application for 
small unattended power sources. If such develop- 
ments could demonstrate the general feasibility of 
the ionic generators, the field of application is 
potentially very wide because of its simplicity, its 
good thermal cycle and small size. 


The author is grateful to the N.R.D.C. for permis- 
sion to mention their interest and work. Opinions and 
conjectures expressed are the author’s own but he is 
very grateful to the Ministry of Power for the oppor- 
tunity to make this study and permission to publish 
this note. 
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Appendix | 

For the convenience of the reader a complete 
alphabetical list of all the symbols used in the 
article is given below. 


A receiving area of detector 

a = absorption coefficient for radiation 

B magnetic flux 

C: constant describing fluctuations in 
photon noise 

E’ constant describing flicker noise 

e electronic charge 

y frequency, modulation frequency 

Af bandwidth 

| d.c. bias current 

I - current corresponding to _ incident 


radiation flux Q, 


Atomic Energy (Geneva, 1958) New York: United 
Nations, 1959 

k Boltzmann’s constant 

Lp diffusion length 

Lpe.pn= electron and hole diffusion length 
Meh electron and hole mobility 

n equilibrium electron density 

Qo mean quantum flux density at surface 
S surface recombination velocity 

ii absolute temperature 

T decay constant of excess carriers 

< carrier transit time between electrodes 
a = bias voltage 

V, = signal voltage 

o = angular frequency = 2z/ 

x inter-electrode distance 

4 specimen thickness, junction depth 

Z impedance 

















MICROWAVE AMPLIFICATION BY 
STIMULATED EMISSION OF RADIATION 


C. R. DITCHFIELD 


Royal Radar Establishment, Malvern 


Stimulated emission, produced by interaction between microwave radiation and a molecular 


system, can result in coherent amplification or oscillation. 


At these frequencies, spontaneous 


emission, which causes random electrical noise, is dominated by stimulated emission. Mech- 

anisms which produce excess noise, being absent, the Maser is an extremely low-noise amplifier. 

The physical processes underlying absorption and radiation, together with the potentialities 
of gaseous and solid state Masers, are described. 


M.A.S.E.R., an acronym for the title of this article, 
was coined by Professor Townes! of Columbia 
University in 1955 to describe the process whereby 
energy, stored in a molecular or atomic system, 
could be used to provide coherent amplification or 
oscillation at microwave frequencies. The subject 
has two main facets: firstly, the provision of ex- 
tremely low noisg¢ amplifiers and very monochro- 
matic oscillators, which has obvious practical 
attractions; secondly, the theoretical interest in a 
process of interaction between radiation and matter 
which is virtually impossible to observe at optical 
frequencies, where the physicist has accumulated 
most of his spectroscopic data. An outline of the 
theory is presented as an aid to the appreciation of 
practical Masers and their future possibilities. 


Absorption Lines 
The quantum theory postulates that an atom (or 
molecule) can exist only in a series of states corres- 
ponding to discrete values of energy. The differ- 
ences between the energy levels, £,, E., E; etc 
decide the frequencies of the radiation which will 
interact with the atom, according to 

hv,, = E, — Ey 

hve, = E; — E, etc 
where h is Planck’s constant, 6-624 x 10-?? erg sec, 
and vz, is the frequency appropriate to a transition 
between levels a and b. If the atom absorbs a quan- 
tum, Av, it jumps from a lower to a higher energy 
level, the reverse process being accompanied by the 
emission of a quantum. 


For a single atom the energy levels are quite 
discrete and sharp lines are obtained. However, in 
aa atomic or molecular conglomeration, particu- 
larly in a solid, interaction occurs between neigh- 
bours, resulting in a broadening of the energy levels 
and a consequent width to the absorption line. 
Instead of a single frequency, a narrow band of 
frequencies can be absorbed (or emitted) and a 
characteristic line shape is obtained, as in Figure /, 


when the magnitude of the absorption is plotted 
against frequency. The details of the shape will 
depend on the broadening mechanism; the Lorent- 
zian line shape, for example, is represented by 
; &(V9) 
1 + {2a(v — v9)T.}* 

In thermal equilibrium the number of atoms in each 
level, i.e., the population, is determined by the 
Boltzmann distribution. Thus for two levels we have 
a population ratio 

Np ~ hey 

=¢ sr 

Na 
where k is Boltzmann’s constant 1-38 x 10-'* ergs 
per degree and 7 is the absolute temperature. For 
microwave frequencies when hy << kT, the exponen- 
tial is often approximated to the first terms of its 
series. 


g(v) 


The probability of a particular atom suffering a 
transition has the same value for either the absorp- 
tion or emission of a quantum. Consequently the 
nett absorption of power for an atomic system is 
proportional to 

W,2(N, — NAV2 

where W,, represents the probability of a transition 
between the levels occurring in unit time. In thermal 
equilibrium the population N, of the lower level is 
the greater and there is a nett absorption of power. 
This absorption disturbs the values of the popula- 
tions, but there is also a tendency for the system to 
relax to its state of thermal equilibrium by some 
system of energy interchange. Hence, for small 
incident powers, the populations do not differ 
significantly from the equilibrium values. However 
if the relaxation time is sufficiently long, a large 
incident power can saturate the transition, that is, 
make N, equal to N,; no further absorption can 
take place and the substance is transparent. This 
process cannot make the population of the upper 
level exceed that of the lower level, because of the 
identity of transition probabilities for upward and 
downward transitions, 


265 











266 C. R. DITCHFIELD 


The observation of absorption lines in the micro- 
wave region has led to the techniques of microwave 
spectroscopy?-‘, which have been used to elucidate 
such problems as that of interatomic bonding. 





Absorption 











Vo 
Frequency 
Figure 1. Lorentzian line shape 





However, no direct use can be made of this absorp- 
tive process to provide amplification of a micro- 
wave signal and, although it is basically capable of 
providing a frequency standard, as discussed by 
L. Essen®, the degree of monochromatism is limited 
by the relatively small population difference be- 
tween the levels. By making the population of the 
upper level exceed that of the lower and thus pro- 
viding stimulated emission of radiation improved 
amplifiers and oscillators can be obtained. 


Stimulated Emission of Radiation 


Visible radiation, for example from a neon discharge, 
is incoherent, being caused by spontaneous emission, 
a random release of quanta as some of the upper 
population fall to a lower energy level. Useful 
amplification demands coherence, whereby each 
emission is stimulated by the electromagnetic wave. 
Einstein has shown that the probability ratio of 
spontaneous to stimulated emission, for thermal 
equilibrium, is given by 
enho/kT — | 

This has a large value at optical frequencies, but is 
small at longer wavelengths; for room temperature 
and a ten centimetre wavelength, the probability 
ratio is about 5 = 10-*, In the microwave region the 





spontaneous emission, which corresponds to electri- 
cal noise, is dominated by the process of stimulated 
emission, which, in a Maser, can provide either very 
low noise amplification or an extremely pure 
spectrum in an oscillator. The problem is to force 
a population excess into the upper energy level. 


The Ammonia Maser 

The molecular beam Maser’, due to C. H. TOWNES 
and his colleagues, was the first to be successful. 
The electric dipole moment of the ammonia mole- 
cule results in a pair of energy levels which provide 
resonant absorption at a wavelength in the region of 
1-25 cm. The required population excess is obtained 
by the action of an inhomogeneous d.c. electric field 
on a stream of gaseous molecules. A physical separa- 
tion occurs, according to the energy state of each 
molecule, so that the higher energy molecules are 
focussed through a hole in a resonant microwave 
cavity. The coincidence of the cavity resonance 
with that corresponding to the energy spacing of 
the ammonia molecules ensures that they are stim- 
ulated to emit radiation, which builds up the micro- 
wave field in the cavity. With sufficient molecules 
present, the molecular power emitted is more than 
enough to overcome the losses in the cavity and 
amplification (or oscillation) is possible. 


In oscillation the ammonia Maser has provided 
the most monochromatic source of radiation ever 
achieved: a spectral purity of a few parts in 10". 
The frequency emitted is characteristic of the 
ammonia molecule and is eminently suitable for a 
frequency standard, the stability being a consider- 
able improvement on that obtained from an absorp- 
tion line. Vibration, chemical and thermal changes 
alter the resonant frequency of the cavity, which in 
turn, detunes the oscillation frequency, though to a 
very much smaller extent. Frequency discrimination 
techniques have been devised to further minimize 
the effect and a long term stability equivalent to 
better than a second a century seems likely. 

As an amplifier the ammonia Maser has the dis- 
advantages of operating at only one frequency and 
of providing useful gain over onlya small bandwidth, 
i.e. a few kilocycles per second, with small output 
power. Some restriction arises inevitably, because 
of the relatively low concentration of molecules, 
thus limiting the energy available for signal amplifi- 
cation. Solids provide a greater concentration of 
active centres and, in particular, use has been made 
of paramagnetic resonance, wherein the micro- 
wave field interchanges energy with the magnetic 
dipole moment of the electron system. Although 
this moment has only about one per cent of the value 
of the electric dipole moment used in the ammonia 
Maser, the higher conceiitration more than offsets 














MICROWAVE AMPLIFICATION 


this and solid state Maser amplifiers have provided 
bandwidths of tens of megacycles per second. More- 
over the resonant frequency of the system can be 
changed by the application of an external magnetic 
field, thus providing a tunable amplifier. 


Paramagnetic Resonance 


The hypothesis of the spinning electron, introduced 
to explain the multiplicity of spectra such as the 
sodium doublet, together with the Pauli exclusion 
principle, can be used to explain the arrangements 
of the orbital electrons round the nucleus of an 
element, leading to the Periodic Table. Associated 
with its spin, s, an electron has a magnetic moment 
uu. The magnetic moment of an ion is the vector 
sum of those of its individual electrons. Our par- 
ticular concern is with the paramagnetic ions of the 
transition groups, titanium to nickel, and the rare 
earths. These possess large magnetic moments be- 
cause of the non-zero spin vector of an incomplete 
intermediate shell of electrons. (Paramagnetism 
can also be obtained from excess electrons associated 
with crystalline lattice defects: those caused by 









































3 
@ = 0° 
2L- 
i- ; 
Magnetic 

a field \, 
ta 2000 4000 

if 

2 

3 

(a) 
3 
6 =60° 
2 = 
Magnetic 

1 field 
> 
e 

0 l i 
us * [700 000 

te 

73 3 

3 

(c) 


Figure 2, Energy levels of Cr+++ 
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neutron bombardment of by 
phosphorus in silicon.) 


impurities, e.z., 


In ions such as trivalent chromium, the total 
orbital angular momentum vanishes and the mag- 
netic moment is that of electron spins, the possible 
values for chromium being + $ and + 4 Bohr mag- 
netons. The chromium ion experiences an internal 
field due to its neighbours in the crystalline lattice. 
This causes the $ values to have an energy different 
to the 4 values, the extent of this crystalline field 
splitting depending on the specific lattice. When 
subjected to an external magnetic field the spin 
energy depends on the sign of the spin, i.e., whether 
its moment has a component parallel or antiparallel 
to the field. This Zeeman splitting thus produces 
four energy levels for chromium. A simple linear 
relation holds with the d.c. field along the crystalline 
axis (Figure 2a), but the picture changes as the angle 
between these two directions is varied (Figures 2b 
to 2d). 


Paramagnetic absorption may occur when the 
frequency of the incident radiation corresponds to 
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in ruby for various angles 0 between the external 
magnetic field and the c axis of the crystal. 
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in kilomegacycles per second to facilitate conversion to transition frequency (E = hv) 
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The observation of absorption lines in the micro- 
wave region has led to the techniques of microwave 
spectroscopy?-‘, which have been used to elucidate 
such problems as that of interatomic bonding. 
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However, no direct use can be made of this absorp- 
tive process to provide amplification of a micro- 
wave signal and, although it is basically capable of 
providing a frequency standard, as discussed by 
L. Essen®, the degree of monochromatism is limited 
by the relatively small population difference be- 
tween the levels. By making the population of the 
upper level exceed that of the lower and thus pro- 
viding stimulated emission of radiation improved 
amplifiers and oscillators can be obtained. 


Stimulated Emission of Radiation 


Visible radiation, for example from a neon discharge, 
is incoherent, being caused by spontaneous emission, 
a random release of quanta as some of the upper 
population fall to a lower energy level. Useful 
amplification demands coherence, whereby each 
emission is stimulated by the electromagnetic wave. 
Einstein has shown that the probability ratio of 
spontaneous to stimulated emission, for thermal 
equilibrium, is given by 
emhv/kT — | 

This has a large value at optical frequencies, but is 
small at longer wavelengths; for room temperature 
and a ten centimetre wavelength, the probability 
ratio is about 5= 10-*, In the microwave region the 





spontaneous emission, which corresponds to electri- 
cal noise, is dominated by the process of stimulated 
emission, which, in a Maser, can provide either very 
low noise amplification or an extremely pure 
spectrum in an oscillator. The problem is to force 
a population excess into the upper energy level. 


The Ammonia Maser 

The molecular beam Maser’, due to C. H. TOWNES 
and his colleagues, was the first to be successful. 
The electric dipole moment of the ammonia mole- 
cule results in a pair of energy levels which provide 
resonant absorption at a wavelength in the region of 
1:25 cm. The required population excess is obtained 
by the action of an inhomogeneous d.c. electric field 
on a stream of gaseous molecules. A physical separa- 
tion occurs, according to the energy state of each 
molecule, so that the higher energy molecules are 
focussed through a hole in a resonant microwave 
cavity. The coincidence of the cavity resonance 
with that corresponding to the energy spacing of 
the ammonia molecules ensures that they are stim- 
ulated to emit radiation, which builds up the micro- 
wave field in the cavity. With sufficient molecules 
present, the molecular power emitted is more than 
enough to overcome the losses in the cavity and 
amplification (or oscillation) is possible. 

In oscillation the ammonia Maser has provided 
the most monochromatic source of radiation ever 
achieved: a spectral purity of a few parts in 10%. 
The frequency emitted is characteristic of the 
ammonia molecule and is eminently suitable for a 
frequency standard, the stability being a consider- 
able improvement on that obtained from an absorp- 
tion line. Vibration, chemical and thermal changes 
alter the resonant frequency of the cavity, which in 
turn, detunes the oscillation frequency, though to a 
very much smaller extent. Frequency discrimination 
techniques have been devised to further minimize 
the effect and a long term stability equivalent to 
better than a second a century seems likely. 

As an amplifier the ammonia Maser has the dis- 
advantages of operating at only one frequency and 
of providing useful gain over onlya small bandwidth, 
i.e. a few kilocycles per second, with sinall output 
power. Some restriction arises inevitably, because 
of the relatively low concentration of molecules, 
thus limiting the energy available for signal amplifi- 
cation. Solids provide a greater concentration of 
active centres and, in particular, use has been made 
of paramagnetic resonance, wherein the micro- 
wave field interchanges energy with the magnetic 
dipole moment of the electron system. Although 
this moment has only about one per cent of the value 
of the electric dipole moment used in the ammonia 
Maser, the higher concentration more than offsets 
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this and solid state Maser amplifiers have provided 
bandwidths of tens of megacycles per second. More- 
over the resonant frequency of the system can be 
changed by the application of an external magnetic 
field, thus providing a tunable amplifier. 


Paramagnetic Resonance 


The hypothesis of the spinning electron, introduced 
to explain the multiplicity of spectra such as the 
sodium doublet, together with the Pauli exclusion 
principle, can be used to explain the arrangements 
of the orbital electrons round the nucleus of an 
element, leading to the Periodic Table. Associated 
with its spin, s, an electron has a magnetic moment 
uu. The magnetic moment of an ion is the vector 
sum of those of its individual electrons. Our par- 
ticular concern is with the paramagnetic ions of the 
transition groups, titanium to nickel, and the rare 
earths. These possess large magnetic moments be- 
cause of the non-zero spin vector of an incomplete 
intermediate shell of electrons. (Paramagnetism 
can also be obtained from excess electrons associated 
with crystalline lattice defects: those caused by 









































3 
6 = 0° 

2b 

1 
= Magnetic 
a field \. 
| 2000 4000 

\- 

2 

3 

(a) 
3 
6 = 60° 
2 = 
Magnetic 

Ir field 
> 
4 

n i 
& . 2000 4000 
i 

1 

2t 

3 

(c) 


Figure 2. Energy levels of Cr+++ 


267 


neutron bombardment 
phosphorus in silicon.) 


ot by impurities, e.g., 


In ions such as trivalent chromium, the total 
orbital angular momentum vanishes and the mag- 
netic moment is that of electron spins, the possible 
values for chromium being + $ and + 4 Bohr mag- 
netons. The chromium ion experiences an internal 
field due to its neighbours in the crystalline lattice. 
This causes the $ values to have an energy different 
to the $ values, the extent of this crystalline field 
splitting depending on the specific lattice. When 
subjected to an external magnetic field the spin 
energy depends on the sign of the spin, j.e., whether 
its moment has a component parallel or antiparallel 
to the field. This Zeeman splitting thus produces 
four energy levels for chromium. A simple linear 
relation holds with the d.c. field along the crystalline 
axis (Figure 2a), but the picture changes as the angle 
between these two directions is varied (Figures 2b 
to 2d). 


Paramagnetic absorption may occur when the 
frequency of the incident radiation corresponds to 
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in ruby for various angles 0 between the external 
magnetic field and the c axis of the crystal. 


Magnetic field in oersteds; energy 


in kilomegacycles per second to facilitate conversion to transition frequency (E = hv) 
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Figure 3. 
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Tne effect onthe signal absorption line of successively higher levels of pump power. 


These photographs were obtained using the R.R.E. ruby Maser working with liquid oxygen as the refrigerant 


the energy gap between two levels; pictorially when 
this frequency coincides with the Larmor precession 
frequency of the electron gyromagnet. Where the 
energy levels are non-linearly related to the magnetic 
field, they cannot be characterized by a single 
quantum number (an ‘admixture of states’) and 
transitions can occur which would be forbidden, 
normally, by the selection rules. This fact is funda- 
mental to the concept of the three level Maser. 


The proximity of a similar neighbour influences a 
paramagnetic ion, the effect being to broaden the 
absorption line as the ions are brought closer to- 
gether. For Maser purposes the active ion is 
‘diluted * by an isomorphous diamagnetic ion, e.z., 
in ruby, an aluminium oxide lattice with one alu- 


minium ion in, say, a thousand replaced by para- 
magnetic chromium. The degree of dilution is a 
compromise between the advantage of having the 
maximum number of active centres per unit volume 
and the disadvantage of having absorption lines so 
broad that the energy levels cannot be treated as 
ideal, discrete and independent states. 


The Two Level Solid State Maser 

In a system with only two levels, i.e., spin + 4, most 
of the spins, N,, will be in the lower energy state 
with their magnetic moment » parallel to the field 
and fewer spins N, will have their moments anti- 
parallel. Various methods of achieving a population 
excess in the upper level can be postulated, the simplest 
to visualize being a reversa! of the d.c. field in a time 
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short compared with the relaxation time. This experi- 
nient waSsperformed by E. M. Purce.t and R. V. 
POUND using nuclear spins, but these have moments 
about a thousand times smaller than those of 
electrons and the stimulated emission available did 
not exceed the circuit losses. 


An easier process to operate is that of adiabatic 
fast passage, wherein the frequency of an oscillator 
is swept rapidly from one side of the resonance value 
to the other. This inverts the individual moments and 
thus has essentially the effect of reversing the popu- 
lation values in the upper and lower states. The 
atomic system is now able to provide stimulated 
emission of radiation, but the gain will decrease as 
the relaxation process tends to restore the popula- 
tions to the values corresponding to thermal equilib- 
rium. Cooling the system has two advantages: it 
increases the relaxation time and it decreases the 
limiting noise temperature of the amplifier. The 
first attempt® to obtain Maser action in the solid 
state used the two levels present in phosphorus 
doped silicon codled to 2°K. The two level Maser 
provides only intermittent operation and much 
more effort has gone into the three level Maser. 


The Three Level Solid State Maser 


In the three level Maser, as first propounded by Pro- 
fessor BLOEMBERGEN’ in 1956, power is applied at 
a frequency corresponding to the transition between 
non-adjacent energy levels, for example, 1 to 3. There 
is a consequential tendency to equalize the popula- 
tions in these two levels, with appropriate changes in 
N,/N. and N,/N,, as the pump power between the 
outer levels is increased. Saturation of the pump 
transition, that is N, = N,, ensures that either 
N3> Nz or N.>N,. The alternative occurring in 
practice, for any particular substance, depends on 
the relative position of the intermediate level and 
on the transition probabilities between the levels. 
These probabilities depend on the strength and 
orientation of the microwave magnetic field and 
also on the transition involved. 


The effect on the signal absorption line (in this 
instance between levels 2 and 3) of successively 
higher pump powers is shown in Figure 3. It can be 
seen that the absorption in the centre of the line 
is decreased until the substance is transparent: 
\. == N, and then, as the ratio N,/N, increases, more 
stimulated emission becomes available, the circuit 
losses are overcome and the system amplifies. 

The process provides continuous operation and 
the noise temperature of the signal transition, 
decided by the ratio N,/N., can be a fraction of 
ambient temperature. Advantage is again obtained 
by cooling the system to decrease noise and to 
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obtain long relaxation times, thus decreasing the 
pump power required, only a few microwatts being 
necessary for ruby at helium temperatures. 


Although the fundamental idea of the three level 
Maser has the beauty of being essentially simple, 
some of the finer details have a complexity which 
has not been elucidated. Nevertheless many three 
level Masers are now working. The main centre of 
effort is the United States, with smaller numbers in 
the United Kingdom, Australia and Holland. Some 
of these achievements will serve to illustrate practical 
considerations of Maser design and to indicate the 
performance to be expected. 


Practical Three Level Maser Systems 


In the first successful three level Masers*-'* the 
interaction between the microwave field and the 
spin system was increased by placing the para- 
magnetic sample in a cavity, which was resonant 
simultaneously to the pump and signal frequencies. 
An external field of appropriate magnitude and 
orientation ensured that the transitions in the 
crystal corresponded to these frequencies. Amplifi- 
cation is possible when the power available from 
stimulated emission exceeds that absorbed in the 
cavity walls; the device oscillates if the stimulated 
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Figure 4. Maser gain and bandwidth as a 
function of input coupling. The experimental 
points were obtained with the first Maser 
working near to 60°K. Better gain-band- 
widths would be obtained at lower temperatures 


emission is sufficient to overcome also the power 
lost through the coupling between the cavity and the 
remainder of the microwave system. Expressed in 
terms of the Q factor, defined as 
O= Bes Stored energy 
power absorbed 
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the expression for stable amplification may be 
written as 


1 1 * l ‘ 1 
Q A \Qm| Qu 
where Q, is the unloaded Q factor of the cavity, 
Q, is the Q factor of the coupling aperture and Q,, 
is the magnetic Q of the sample. Q,, is negative in 
sign for an emissive sample. 


To achieve good performance |Q,,| should be as 
low as possible, as this determines the bandwidth 
and gain which can be obtained. For a reflection 
cavity Maser, in which the incident signal is ampli- 
fied and reflected back along the input line in the 
reverse direction, the formulae for bandwidth B, 
and power gain G, are 


Ba (3; 4a - o.)) 


2 


1 l l 
a a, *(1ea- 3) 


l l 1 
Q, (ioe ac} 


Therefore the quantity 


ic l 1 
s|ve-1|=>| 19-9, 


is a constant for a particular system. Thus the 
bandwidth may be increased, with a consequent 
decrease in gain, by increasing the loading on the 
cavity, that is decreasing Q,. This is shown graph- 
ically in Figure 4. The system described has the 
disadvantage that it requires a circulator to distin- 
guish between the incident and amplified signals. A 
transmission cavity Maser, for which essentially 
similar expressions can be derived, requires only an 
isolator at its output to prevent noise from subse- 
quent components in the system reaching the Maser 
and being amplified. However its gain-bandwidth 
product for the same value of |Q,,| is somewhat 
worse than the reflection cavity Maser, which has 
been preferred in experiments to date. 


Only a few months after the publication of 
Bloembergen’s proposal for a three level Maser, 
workers at the Bell Telephone Laboratories, an- 
nounced* oscillations at 9-06 kMc/s using a pump 
power of 17-52 kMc/s. The salt used was gadolinium 
ethyl sulphate at 0-5 per cent dilution in a lattice of 
lanthanum ethyl sulphate, and was contained in a 
resonant cavity immersed in liquid helium at 1-:2°K, 
with a magnetic field of 2850 oersteds. A rectangu- 
lar cavity was used resonant in an Hp», mode at 
17-5 kMc/s, and a strip line mode at 9-06 kMc/s 
with the salt placed as shown in Figure 5. 
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On a simple theory of the electron spin gyro, a 
transition between adjacent levels has the maximum 
probability when the r.f. field is perpendicular to the 
d.c. field, as it can couple most effectively with the 





Figure 5. A Maser cavity (by courtesy of 

The Editor of Bell Laboratories Record) 
component of spin moment precessing at the 
Larmor frequency. In general, with the admixture 
of quantum states, the transition probabilities for 
signal and pump frequencies are a maximum at 
particular orientations of the respective r.f. fields 
with relation to the d.c. field. These optima are 
calculable from the matrix elements of the wave 
functions corresponding to the states of the para- 
magnetic ion. A rectangular cavity has the advant- 
age that the direction of the r.f. field can be 
substantially the same over quite a large sample. 
Therefore spins in all parts of the crystal will react 
most efficiently with the pump and signal r.f. fields, 
particularly if the transition probability requirements 
are for parallel or orthogonal fields, which can be 
conveniently arranged in such a cavity. 

In the Bell Telephone Laboratories cavity the r.f. 
fields would have been optimum when orthogonal, 
whereas they were, in fact, parallel. Consequently 
the d.c. field was arranged to be at 45° to the rf. 
fields, so that the signal field had a component 
perpendicular and the pump field a component 
parallel to the external magnetic field. In a cavity of 
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circular cross section the direction of the field will 
change over the sample and only a component of it 
‘s really effective. Thus, two factors which affect 
Qm| are the orientation of the pump field, which 
affects the extent to which this transition may be 
caturated and the orientation of the signal field which 
determines the ease with which stimulated emission 
is obtained from the spins. 


A filling factor is defined as the proportion of the 
energy in the cavity which would be stored in the 
volume occupied by the sample, that is 

SH?dv 

fH?dv 
where v and V refer to the sample and the cavity 
volumes. The pump and signal fields will vary in 
intensity throughout the cavity to a different extent 
and it would appear that the advantage of a large 
filling factor is offset if part of the crystal contains a 
node at the pump frequency, as, in this region, the 
populations would not be disturbed. In practice 
completely filled»cavities have been used, a possible 
explanation being that phonons can be responsible 
for energy interchange between different parts of the 
crystal, thus resulting in an adequate population 
ratio over the whole volume. 


The final contributions to |Q,,| arise from the 
available spin population per unit volume per unit 
line width, effectively N7., and their distribution to 
give a large ratio between the populations in the 
states apertaining to the signal transition. Because 
the state populations decrease exponentially with the 
energy difference, i.¢., in proportion to the frequency 
of transition, there is advantage in using a high value 
of pump frequency, as, other things being equal, 
this will produce a large ratio of the signal level 
populations. 


A very interesting variation on this is the method 
of ‘ push-pull pumping’, which can be applied 
particularly when two suitable pump frequency 
transitions occur at the same frequency, for a given 
magnetic field. For example, at a particular angle 
of 54° 44’ between the d.c. field and the crystalline 
axis of ruby, the energy levels are symmetrical, as 
in Figure 6. Thus pump power applied between 
levels 1 and 3 will also act between levels 2 and 4, 
pushing spins into level 3 and pulling them out of 
level 2. The result is a value of |Q,,| which is 
several times lower than can be obtained by simple 
pumping between one pair of levels. With other 
ions such as Fe** in sapphire, degeneracy can occur 
between more than two transitions, thus giving 
further scope to this idea. 


1 


So far it has been assumed that the width of the 
ibsorption line, expressed as a frequency, is much 


greater than the amplifying bandwidth but with 
more recent developments in practical Maser sys- 
tems this is not so. Therefore, in addition to Q,,| 
a factor Q, apertaining to the absorption line and 
defined as 27T,, must be taken into account. If 
|Qm| is very low, Q, is the limiting factor in band- 
width performance and it could be advantageous to 
broaden the absorption line, that is to decrease T,, 
by use of an inhomogeneous d.c. field, thus adjusting 
Q, to be approximately the same value as |Q,,,. 


Cavity Masers have been constructed which are 
tunable’, or can operate at relatively high tempera- 
tures, or can provide gain bandwidth products” of 
nearly 100 Mc/s, that means, a bandwidth of several 
megacycles per second with a gain of 20 dB. Further 
research, particularly into paramagnetic materials, 
might lead to all these features being combined in 
one device, with perhaps only a very small external 
magnet, particularly if alarge splitting can be obtained 
from the crystalline field. 
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Figure 6. Energy levels of Cr+++ in ruby. 
Magnetic field at 54° 44’ to c-axis. (Energy 
expressed in kMc/s) 


However, another type of microwave circuit, the 
travelling wave structure, has many advantages over 
a cavity e.g. it has good gain stability, tunability, 
its amplifying bandwidth is less limited by Q, and 
the unidirectional property, provided by an external 
device in a cavity system, can be built into the Maser 
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itself. Such a Maser has been built’ by Bell Tele- 
phone Laboratories (Figure 7). The wave velocity 
is decreased by a comb-like structure which effect- 
ively increases the interaction between the wave and 
the crystals. The circular polarization of the wave, 
which has opposite sense either side of the comb, is 
cleverly utilized to give amplification in one direc- 
tion, when the wave interacts with the Larmor 
precession of ruby with a low chromium concen- 
tration, but attenuation for a wave in the other 


Emissive | Absorptive 











Figure 7. Diagram of travelling wave struc- 
ture in a ruby Maser used at Bell Telephone 
Laboratories 


direction, which interacts with heavily doped ruby, 
in which the populations are not brought into the 
emissive condition by the pump power. A gain of 
about 25 dB over a bandwidth of 25 Mc/s has been 
obtained. 


Noise in Masers 


It has been customary to judge receivers by their 
noise factor. This is an adequate criterion for 
relatively poor receivers, e.g., of noise factor six 
times or more, where the output noise is little 
affected by the background temperature against 
which a signal is to be detected. However, with low 
noise amplifiers, such as Masers, the reverse is true 
and a much better criterion is the effective noise 
temperature contribution of the Maser amplifier. 
It can be shown that 
a 


4 
Tey =- re: ¢  « [Tul —* 
Q\2ni| 6. on a, 


for a reflection cavity, where 7, is the temperature 
of the cavity and T7,, is defined in terms of the 
emissive populations of the signal transition energy 
states. A large ratio decreases |7,,,| as well as |Q,,,\, 
thus emphasizing again the advantage of such 
schemes as push-pull pumping. These can result in 
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|Tm| being below ambient temperature. The 
sensitivity limit of a three level Maser can thus be 
below the physical temperature of the device. 


In practice the noise performance of a Maser is 
so good that it is extremely difficult to measure accur- 
ately, being about a thousand times lower than that 
of aconventional receiver. However an upper limit 
of 10° to 20°K has been set by A. L. MCWHorTER, 
J. W. MEYER and P. D. Strum!® of Lincoln Labora- 
tories on a cavity Maser using potassium chromi- 
cyanide and by R. W. DeGrasse, E. O. ScHULz- 
DuBois and H. E. D. Scovit, of Bell Telephone 
Laboratories, on the travelling wave ruby Maser. 


Thus there is little doubt that a system using a 
Maser will be limited by noise contributions from 
other parts of the equipment and work will be 
necessary to improve aerials and transmission 
circuits from this point of view. 


Application of Masers 


To date Masers have been developed mainly for 
their own interesting properties, with obvious appli- 
cation in the future to low noise amplifiers used in 
scientific investigations such as microwave spectro- 
metry and radio astronomy. In the latter field 
the Naval Research Laboratories, Washington, 
in cooperation with Columbia University, have 
used a Maser as a pre-amplifier for a _ radio 
telescope, with great advantage over a con- 
ventional receiver. Other applications lie in the 
realms of radar and long distance communication, 
perhaps the most exciting being those techniques 
concerned with global communications via satellites 
and, eventually, inter-planetary travel. Finally, 
although all present systems involve use of a pump 
frequency higher than that at which stimulated 
emission is obtained, several possibilities exist 
whereby coherent oscillators and amplifiers may be 
obtained at submillimetre wavelengths, which region 
is at present undeveloped. 


The extent to which Masers are used will depend 
on the extent to which their low noise properties can 
be utilized effectively in asystem and on the degree of 
complexity in the Maser itself. At present the pro- 
vision of useful gain over large bandwidths entails 
very low temperatures but, in a multi-level Maser, 
this in not fundamentally necessary. A large popu- 
lation ratio for an emissive signal transition can be 
obtained by pumping several pairs of levels and by a 
favourable adjustment of the relaxation times 
between specific energy states. The latter technique 
would also improve the power handling capacity of 
a Maser, which, although completely adequate for 
passive receivers, has so far limited application to 
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radar. In gadolinium ethyl! sulphate one transition 
time was changed by an order of magnitude® on 
introducing a small number of cerium ions into the 
lattice. 


Although data exist on the paramagnetic proper- 
ties of many materials, the precise mechanism of the 
spin-spin and spin-lattice relaxation processes needs 
further elucidation. This could lead to the synthesis 
of a material with desired properties, which, ‘ tail- 
ored’ for Maser applications, might transcend 
present crystals to such an extent that low noise, 
adequate gain and wide bandwidth will be obtain- 
able from a relatively simple Maser working at 
easily achievable temperatures, such as 77°K (liquid 
nitrogen), with only a very small external magnetic 
field for purposes of tuning. 


In conclusion it is apparent that Maser techniques 
have developed rapidly in the few years since their 
conception and it would appear that the ultimate 
possibilities still remain to be exploited and divi- 
dends should continue to accrue from further 
research. . 

This article is published by kind permission of the 
Director of the Royal Radar Establishment and the 
Ministry of Supply. 
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Book Reviews 


Preparation of Single Crystals 

W. D. Lawson and S. NIELSON 
(vii+255 pp; 8} in. by 5% in.) 

Butterworths Scientific Publications: 
York: Academic. 45s; $8.80 


Tuis book is one of a projected series of semiconductor 
monographs and is concerned with the growth of semi- 
conducting crystals. The general recognition of the 
importance of crystal purity and perfection in the 
experimental study and application of these solids, and 
the consequent interest in the improvement of techniques 
of preparation, make this a subject of great current 
interest. 

An initial elementary outline of some electrical 
properties of these materials is followed by a survey of 
methods of crystal growth and a description of equip- 
ment. The authors then discuss means of purification 
and analysis, and procedures for the preparation of 
compounds. 

An account of equilibrium and non-equilibrium 
defects leads to a discussion of the mechanisms of 
growth, together with some further techniques. A final 
chapter, with a title open to misconstruction, is mainly 
concerned with the interaction of impurities with the 
lattice. 


London: New 


See also 


The book has an air of informality, makes very few 
demands on the reader in the form of knowledge of 
semiconducting systems, and provides a clear and 
valuable account of the techniques. The authors 
sometimes go to great lengths to ensure that practical 
and constructional details are given. Many of these are 
very valuable but it does not seem possible to justify 
the inclusion of some of the material. For example, 
in the chapter on apparatus, the elementary theory 
of the McLeod gauge and diagrams illustrating the 
principles of operation of rotary backing pumps are out 
of place. 


More generally, the selection and emphasis of the 
contents may in some way restrict the value of the book. 
The primary aim is clearly to provide a practical account 
of techniques, particularly for the newcomer to the 
subject. This aim is very well satisfied. The authors do 
subsequently give a simple but satisfactory discussion of 
growth mechanisms, but do not analyse their objectives. 
Nor do they discuss the design of procedures to produce 
given physical properties. In the field of semiconductors 
the conflicts between different practical objectives or of 
criteria of perfection, would surely have been of interest 
to the reader, however practical his approach, and one is 
sorry that these have not been discussed, P. C. BANBURY 
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THE THERMODYNAMIC 
BEHAVIOUR OF SOLIDS 


M. A. JASWON 


Department of Mathematics, Imperial College of Science and Technology, London 


In addition to the four articles on crystal physics by Dr Jaswon published * 
during 1958, a fifth has now been added. This provides an outline of the main 
principles of thermodynamics, particularly in relation to the light they throw 


on the behaviour of solids. 


(The five articles on crystal physics will shortly be 


reprinted in the form of a paper-covered booklet.) 


A THERMODYNAMIC system can exist in several 
distinct kinds of state relative to its environment, 
but they fall into two main classes, equilibrium states 
and non-equilibrium states. The first possibility 
means that the system has no tendency to change 
with time in any way, though to check this rigor- 
ously would require an exhaustive series of experi- 
ments. For instance, a metal rod might appear to 
be in equilibrium at a particular temperature and 
pressure, as judged by the fact that it persisted 
indefinitely as such. However, if it had previously 
been cold-worked, or quenched, or subjected to 
radiation damage, its microcrystalline structure 
would certainly be changing for some time after- 
wards. The change of internal structure can some- 
times be followed by x-ray analysis, as shown by 
Figure 11 (page 19) referring to the self-recovery 
at room temperature of cold-worked silver filings. 
This process of recovery essentially involves the 
elimination of inhomogeneous microstress fields, and 
therefore the introduction of increasing regularity 
throughout the lattice, until a limiting degree of 
order has been achieved. Cold-worked metal thus 
has the status of a non-equilibrium system relaxing 
towards equilibrium at a slow but steady rate, over 
an appreciable time interval. Similar considerations 
apply to the elimination of quenched-in lattice defects 
or quenched-in disorder and, in fact, to any process 
involving the nucleation of an equilibrium system 
at low temperatures. Generally speaking, the reverse 
process of high temperature nucleation takes place 
much more rapidly. For instance a block of ice 
suddenly heated beyond 0°C will immediately trans- 
form into water, so that the non-equilibrium phase 
—ice above the melting point—has only a fleeting 
existence. 


Perhaps the most important practical example of 
a non-equilibrium system is provided by glass, 
which possesses a non-crystalline structure of some- 
what ill defined character. According to the recent 
geometrical considerations of J. D. BERNAL’, a 
* Research, Lond. (1958) 11, 26, 108, 227, 315 


large number of equal spheres can be close packed 
in two—and only two—fundamentally different 
ways. These conform to long range and short range 
order respectively. No intermediate degree of order 
is admissible. If so, the existence of two condensed 
phases of matter becomes understandable and we 
may also infer that the atomic structure of glass 
must be closely allied to the structure of liquids. 
Yet, paradoxically, many properties of glass are 
much more reminiscent of crystalline solids than of 
liquids. For instance, Poisson’s ratio has the values 
0-245, 0-270 and 0-5 for glass, steel and liquids 
respectively. A more direct indication is given by 
the specific heat curve, sketched in Figure | for 
glycerol. As the temperature of the melt, or super- 
cooled liquid, is reduced through a small interval 
about the transition temperature T,, the specific heat 
sharply declines to a value not very different from 
that of the crystalline phase. 


This behaviour at first sight admits to two inter- 
pretations, of almost opposite significance. It could 
indicate an equilibrium order—disorder transition 
of some internal degree of freedom, ¢.g., molecular 
rotation, analogous to the various lambda transit- 
ions in solids. This possibility has been considered 
unlikely by Sir FRANCIS SIMON and his school?. In 
any case such a view implies the almost untenable 
idea of the glassy phase as a possible equilibrium 
phase of condensed matter that competes with the 
crystalline phase down to the lowest temperatures. 
On the second interpretation, now generally accept- 
ed, the fall in specific heat indicates the cooperative 
freezing-in of whatever order has been attained by 
the system at 7,. Attainment of further order, in 
the direction of increasing crystallinity, is almost 
completely inhibited below 7,, so that the system 
persists in its non-equilibrium state. More tech- 
nically expressed, the configurational entropy of the 
system is not progressively eliminated on lowering 
the temperature below 7,, but remains locked up 
in the system even at 0°K. This behaviour stands in 
marked contrast to that of equilibrium systems 
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where, by virtue of the Nernst heat theorem (third 
iaw of thermodynamics) the entropy necessarily 
vanishes at 0°K. 


A two dimensional model of a crystal is drawn in 
Figure 2. This depicts the motif unit as a pair of 
atoms linked by a strong, localized bond of covalent 
or ionic character. Neighbouring units, whose 
position and orientation must conform to the trans- 
lational periodicity of the crystal lattice, are linked 
by secondary, weaker bonds of the van der Waals 
type. The main difference between a crystal and a 
glass springs from the fact that neighbouring units 
in the latter do not conform to any lattice symmetry, 
but remain fixed at substantially the configuration 
appropriate to a supercooled liquid at temperature 
T,. This deviation from regularity increases the 
potential energy of the structure, but the increase is 
energetically tolerable since it affects only the second- 
ary bonds. Evidently, however, only molecular 
units of sufficient complexity, capable of interacting 
by a wide variety of weak coupling mechanisms, have 
a chance of forming a disordered array of the 
requisite stability. 


No equation of state, such as that connecting the 
temperature, pressure and volume of a crystal, 
exists for a glass. In the first place, the internal 
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Figure 1. Temperature variation of specific 
heat of glycerol (after Simon and Lange, cf. 
Reference 2) indicating freezing temperature 
T, and glass transition temperature T,— 
the dashed lines indicate expected specific 
heat if glassy phase did not intervene at T, 
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structure of the latter, with derivative physical 
properties, is not uniquely defined by the tempera- 
ture and pressure, but depends on the whole pre- 
vious thermal and mechanical history of the speci- 
men. In the second place, relaxation effects can 
never be entirely absent in a non-equilibrium system, 
with the result that its physical properties alter over 
long intervals of time. On heating a glass, the 
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thermal energy of the molecules soon becomes 
comparable with the interaction energy of the 
secondary bonding forces, so that they acquire 
greater potentialities for relative displacement. 
Eventually the molecules become sufficiently mobile 
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Figure 2. Two dimensional model of 
(a) crystal and (b) glass 


















































for the substance to behave identifiably as a liquid, 
though the precise point at which it ceases to be a 
solid cannot be sharply determined. It may be 
remarked that no generally accepted, quantitative 
theory either of glass or liquids has yet been form- 
ulated. 


The Entropy Principle 
When a system undergoes physico-chemical change, 
its internal energy usually also changes. Part of this 
change is balanced by a flow of heat between the 
system and its surroundings; the remaining part is 
balanced by mechanical interaction. This is the 
first law of thermodynamics. In the usual symbolism, 
dU=dQ-+ dW, ask 
where all the energy increments are formally positive 
and acquired by the system. If the volume of the 
system changes by dV under an external pressure P, 
then dW equals (— PdV). 


Considerations of energy balance do not by them- 
selves suffice to establish the direction of physico- 
chemical change. This role in thermodynamics is 
assigned to the entropy function S. A direct evalu- 
ation of this function lies beyond the scope of classical 
thermodynamics, which is limited to the theorem 
that entropy flow dQ/T accompanies a heat flow 
dQ. Let us now suppose that S may, by some means, 
be computed for two neighbouring states of a 
system: the increment dS then accounts for the net 
change in entropy of the system arising from the 
various interacting factors—configurational order- 
ing or disordering of the atoms, volume and shape 
deformations, heat flow—that make up the physico- 
chemical process in question. Heat flow entropy is 


certainly included in dS, but does not generally 
comprise the whole of dS, for entropy may be 
created or destroyed within the system quite apart 
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from that passing through the boundaries. On the 
other hand, the corresponding entropy change in 
the surroundings can be nothing other than (— dQ/T). 
We now formulate the second law of thermo- 
dynamics as 

dS — dQ/T = 0, awe cee 
which means that the net change in entropy of 
system plus surroundings, considered as a single 
physical unit, must be either positive or zero. 


The equality sign in (2) defines reversible proces- 
ses. In effect, entropy flow here accounts for the 
whole of dS, so that reversing the flow reverses the 
sign of dS and hence the direction of change. 
Alternatively expressed, reversible processes merely 
involve a redistribution of existing entropy between 
the system and surroundings, though with no pre- 
ferred flow in any particular direction: the system 
thus remains in equilibrium with its surroundings at 
all stages of the operative process. 


Reversible changes do not occur in nature but 
are important theoretically, sometimes enabling dS 
to be computed for a given change of state. Consider 
for example the isothermal expansion of a perfect 
gas from V to V+-dV. We suppose heat is absorbed 
by the gas and entirely converted into work done 
by the gas in expanding against the external pressure. 
The first law of thermodynamics thus reads 

dU=dQ+dw=0, 

with dW = — PdV = (— RT/V). dV, 
by virtue of Boyle’s equation of state. For reversible 
changes, we have 

dS — dQ/T = 0, i 
whence dS = dQ/T = — dW/T=(R/V).dV ..(4) 
From (4), the total change in entropy corresponding 
to a finite change in volume at constant temperature 
T may readily be calculated. This simple example 
illustrates an important technique of classical 
thermodynamics: dS is evaluated indirectly by 
envisaging a reversible process which connects the 
relevant final and initial states of the system. 


The inequality sign in the second law of thermo- 
dynamics (2), regulates the direction of what are 
termed natural or irreversible changes*. As will 
appear later, however, it yields little information 
about the rate of relaxation towards equilibrium; 
nor can it predict the particular route of relaxation 
out of several possibilities which may exist. Perhaps 
the simplest example of an irreversible process 
concerns the isothermal expansion of a perfect gas 
thermally insulated from its surroundings. We 
suppose the volume of the container to be suddenly 
expanded from V to V+dV, the external pressure 
being reduced in accordance with Boyle’s law so as 
to maintain 7 constant. The gas molecules will at 
once irreversibly occupy the extra available volume, 
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until they are spread uniformly throughout th« 
whole volume. Now bearing in mind that the samc 
change of state may be achieved reversibly by the 
method of the last paragraph, we see that dS for 
the present process is provided by (4). On the other 











Figure 3. Particle at potential height ¢, 
will tend to slide down well to equilibrium 
configuration at $» 


hand, since no heat flows in or out of the system, 
Equation 3 is simply replaced by the inequality dS 
is greater than zero, where dS happens to be the 
net amount of entropy created by virtue of the 
natural expansion process. 


Other examples, such as the melting of ice, can 
only be discussed qualitatively. Here the entropy 
function S increases owing to the replacement of 
long range order, characteristic of a crystal, by short 
range order characteristic of a liquid. At the same 
time U increases owing to the rupture and bending 
of interatomic bonds: if so, neglecting dW, we have 

dU= dQ >0, 
implying an entropy flow dQ/T into the system. 
On balance, (dS—dQ/T) exceeds zero above 0°C. 


It is important to understand that S may decrease 
during an irreversible process, provided sufficient 
entropy flows into the surroundings to provide a 
compensating increase. Consider the self-recovery 
of a cold-worked metal. Here S decreases owing to 
the restoration of lattice perfection; at the same 
time, the relief of strain energy implies that 

dU=dQ < 0, eo) 
showing that entropy flows into the surroundings. 
We thus have dS negative and ( — dQ/T) positive, 
with (dS — dQ/T) positive on balance. 

As a second example, crystallizations from the 
vapour must always involve an enormous decrease 
in entropy. Part of the decrease arises from the 


volume contraction on solidification and this is 
balanced by an entropy flow into the surroundings 
by much the same mechanism as in the reversible 
isothermal compression of a perfect gas. The second 
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oart of the decrease arises from the build up of 
-onfigurational order and this is compensated by the 
formation of interatomic bonds, resulting in a 
relation of the form (5). Accordingly, whilst dS is 
large and negative, (—dQ/T) is large and positive, 
with (dS—d@Q/T) positive on balance. It has often 
been remarked that the growth of biological organ- 
ization implies a decrease in entropy: all the avail- 
able evidence indicates that an appropriate gain in 
entropy occurs elsewhere*®. If so, it must be 
accepted that the existence of living organisms does 
not contravene thermodynamic principles. 


Mathematical Formulation 


To make further progress, we write 
dS — dQ/T = —dd¢/T > 0, ee) 
where ¢ stands for a thermodynamic function 
having the dimensions of énergy. The introduction 
of this function brings two advantages: it replaces 
increase in entropy by the more familiar concept of 
decrease in energy, and the quantity | d¢é| in some 
sense measures the driving energy for the reaction 
at any stage, since dé equals zero at equilibrium. 
A valuable analogy for non-equilibrium systems is 
provided by a particle placed at the side of a poten- 
tial well (Figure 3). Just as the particle tends to 
slide down the well, thereby diminishing its potential 
energy until equilibrium is achieved, so a non- 
equilibrium system tends to pass through a sequence 
of non-equilibrium states of gradually diminishing ¢. 
If a system starts at ¢, and ends up in equilibrium at 
do, the total driving energy for the transition is 
pr a. po. 
An explicit expression for ¢ depends on the 
circumstances. In the case of a single homogeneous 
phase, its thermodynamic behaviour is completely 
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Figure 4. Typical dependence of free energy 
@ and internal energy U on temperature 
(schematic) 
controlled by a pair of independent variables 7,P 
or 7,V. With the former pair, d becomes identified 
as the Gibbs free energy G, or thermodynamic 
potential, 
G= U—TS+ PV ee 
of the system. Here U, S, V are all, in principle, 
functions of temperature and pressure, and 


therefore vary if the latter vary. Accordingly, for 
given increments d7 and dP, straightforward 
differentiations and collecting of terms gives 
dG = dU + PdV — TdS — SdT+ VdP ....(8) 
= (dQ — TdS)— SdT+ VdP, ....(9) 
on bearing in mind that dU+PdV equals dQ. 


Coupling Equations 2 and 9, we draw the following 
conclusions. Firstly that for equilibrium at constant 
temperature and pressure, 

dG = dQ — TdS= 0 
and secondly for irreversible processes at constant 
temperature and pressure 

dG = dQ — TdS < 0. 
Strictly speaking, this inequality refers to a non- 
equilibrium function G which has not yet been 
explicitly defined. It is supposed, however, that the 
expression (7) still holds, but that U, S and V depend 
not only on temperature and pressure but on 
additional parameters relating to the processes 
taking place. For reversible processes involving 
varying temperature and pressure, the equation for 
dG becomes 

dG = — SdT+ VdP. 
As a corollary, the equilibrium V corresponding to 
a prescribed temperature 7 and pressure P is 
determined by the ‘ equation of state’ 


V = = -...10) 
oP 


For irreversible processes involving a_ varying 
temperature and pressure, 
dG < — Sd7T + VdP. 


If the independent variables are T and V we 
identify ¢@ as the Helmholz free energy, or free 
energy 

F=U-TS .. (11) 
of the system. For given increments d7 and dV, 
dF = dU — TdS — SdT 
dU + PdV —TdS — SdT — PdV 
(dQ — TdS) — SdT — PdV. ae 


From Equations 2 and /2, it may be concluded that 
for equilibrium at constant temperature and volume 
dF=dQ—TdS=0 
and, in the case of irreversible processes at constant 

temperature and volume, 
dF=dQ-— TdS< 0. 


As before, this inequality refers to a non-equili- 
brium F formally defined by (//), but where U and 
S depend not only on temperature and volume but 
on additional parameters relating to the processes 
taking place. For reversible processes involving 
varying temperature and volume, 

dF SdT — PdV. 


As a corollary, the equilibrium P corresponding 
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to a prescribed temperature and volume is deter- 
mined by the equation of state 
a at 
dV) + 
and this equation must, of course, be equivalent to 
Equation 10. Finally for irreversible processes at 
varying temperature and volume 
dF < — SdT— PdV. 
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Figure 5. Free energy and internal energy 

temperature dependence of face-centred cubic 

lithium (¢,,U,) and body-centred cubic 

lithium (¢., U2). The two free energy curves 

cross at a transition temperature T, in the 

neighbourhood of 80°K. Note the jump from 
U, to U, at T, (schematic) 


Phase Transformations 


For most changes involving only solids, the 
product PV remains practically constant in com- 
parison with U and TS, so that no effective distinc- 
tion arises between F and G. Alternatively expressed, 
dW is smal! by comparison with dQ, so that dU 
approximates to dQ; if so, the relation (6) becomes 
dS — dU/T dd/T > 0, 

whence @= U-—TS «sot 
for either temperature and pressure or temperature 
and volume as the independent variables. Now 
bearing in mind that neither U nor S depend very 
critically on temperature, we draw the following 
qualitative conclusions from Equation 14: 


(i) @ approximates to U as T approaches zero, 
since TS then becomes negligible compared with 
U and ¢ approximates to (— TS) as T gets large, 
since TS then dominates U. As a corollary, the 
dependence of ¢ on T must be as shown in Figure 4. 

(ii) Mathematically, d is minimized at a given 
temperature and volume, or temperature and pres- 
sure, by making U as small as possible and S as 
large as possible. From the physical point of view, 
however, these are seen to be mutually incompatible 
requirements: a small U implies mechanical stabil- 
ity, whereas a large S implies configurational 


disorder. In practice, the difficulty is overcome by 
the appearance of structures which emphasize the 
dominant factor at the temperature concerned. 
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(iii) Taking (i) and (ii) together, we see that 
mechanical stability is preferred at low temperatures, 
hence the appearance of solids. At high temper- 
atures, configurational disorder is preferred, 
hence the existence of gases. At intermediate 
temperatures, the two factors are more evenly 
balanced, resulting in the existence of the liquid 
phase as a genuine compromise between them. 


If a substance has a choice of two different crystal- 
line modifications, of slightly differing free energies 
(allotropy, polymorphism), the more close-packed 
structure will be preferred at low temperatures and 
the more open structure as the temperature increases. 
These qualitative considerations led C. ZENER to 
postulate the existence of face-centred cubic lithium 
and sodium at sufficiently low temperatures, a fact 
subsequently discovered by C. S. BARRETT and co- 
workers. Denoting the face-centred cubic and 
body-centred cubic free energy curves by ¢, and 
¢, respectively, their relative behaviour is indicated 
in Figure 5: both have the same general form, but 
cross at a transition temperature 7, at which ¢, 
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Figure 6. Schematic P, T phase diagram. 

D marks the triple point and E the critical 

point—the significance of A, B, C is 
explained in the text 
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Figure 7. Temperature dependence of 

internal energy and specific heat in neigh- 

bourhood of critical temperature T, for 
order—disorder transition (schematic) 
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Figure 8. (a) Model of perfect order: successive lattice rows are alternatively black and 
white; (b) partial order: successive lattice rows are alternatively mainly black and mainly 
white; (c) Disorder: random distribution of black and white in each row 


equals dn, i.e., 

U, — T.S, = U, — T,S:, <iea ee 
and L=U, — U,= TAS, — S)>0,_ ...-(6) 
where L stands for the latent heat of transformation. 
The relation (/6) governs the way in which U and 
S jump at the transition temperature, it being noted 
that both become ‘greater for the high temperature 
phase. 


Quantitative calculations of ¢, and ¢, for sodium, 
undertaken very recently by Dr M. A. E. NutKins® 
on the basis of Debye’s theory of solids, have 
provided a theoretical check on Zener’s ideas and a 
transition temperature of the right order of mag- 
nitude. At temperatures in the neighbourhood of 
the transition temperature 7,, we have the useful 
approximation 

gd: — $; = (U, — U,) — T(S, — S,) = L — TL/T, 

Li — 7/T-), assole 
for the free energy difference between the two phases. 
Similar considerations apply to other polymorphic 
transitions, and to the solid—liquid transition. Solid 
gas and liquid-gas transformations involve such 
enormous volume changes that the Gibbs free energy 
G differs significantly from the Helmholz free 
energy F, thus making it necessary to specify 
whether temperature and volume or temperature 
and pressure are the controlling variables; other- 
wise, the preceding general ideas still hold good. 


It is usual to assemble information about phase 
relations in a pressure-temperature diagram as 
typified in Figure 6. The field of representative 
points divides into a number of regions, each 
characterized by a specific equilibrium phase; two 
chases co-exist in equilibrium along the boundary 
between them, and three phases co-exist at the triple 
point. The liquid and gas phases are conventionally 
distinct only up to the critical point, at which their 
pecific volumes become equalized. Recent work, 
however, suggests that the boundary should be 
ontinued as indicated by the dashed line, the 


transition across the latter being of a type discussed 
in the next section’. 

As an exercise in the use of the phase diagram, 
let us examine what happens to a quantity of liquid 
placed in an evacuated rigid container held at a 
temperature 7,. Not being subjected to any pressure 
initially, the representative point A lies initially on 
the 7-axis at a height 7,: from the diagram we see 
that gas is the equilibrium phase at A, and hence 
deduce that the liquid irreversibly evaporates. As 
evaporation proceeds, P necessarily increases, and 
the representative point therefore moves parallel to 
the P-axis along the dotted line as shown. There 
are now two possibilities. The liquid has been 
completely vaporized at a point B within the gas 
phase field. If so, the system at B has the status of a 
gas of fixed volume (that of the container) currently 
at a temperature 7,. Alternatively, the liquid has 
not been completely vaporized before the phase 
boundary is reached at C. If so, liquid and vapour 
co-exist in equilibrium at C, and transformation has 
no tendency to proceed further. On raising the 
temperature beyond 7,, the representative point 
moves along the phase boundary in the direction of 
increasing pressure, indicating further evaporation. 


Second Order Transitions 

The transitions so far discussed are termed first- 
order since they involve a discontinuity in the 
internal energy U. There exists a class of transform- 
ations which involve no discontinuity in U, only in 
the specific heat dU/d7, and which are hence said 
to be second order. Sometimes they are termed 
lambda-transitions, in view of the characteristic 
shape of the specific heat curve at the critical 
temperature 7, (Figure 7). Perhaps the simplest 
example of such a transition is the order-disorder 
change in binary alloys*. Below 7,, the two kinds 
of atom are arranged perfectly regularly on a 
common lattice; as T approaches 7,, the regularity 
begins to disappear, at first gradually and then 
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catastrophically, until just beyond 7, the distribut- 
tion has become completely random (Figure 8). 


To describe this situation within the framework 
of classical thermodynamics, we introduce an order 
parameter ¢, having the same status as P, V or T: 
perfect order is defined by £ being unity, perfect 
disorder by ¢ being zero, and states of partial order 
by intermediate values of £. As the temperature 
varies in the critical range, £ varies continuously 
from one of its extreme values to the other, as 
sketched in Figure 9. It is convenient here to think 
of each intermediate structure, for example that 
characterized by ¢,, as existing hypothetically over 
the whole temperature range, with a free energy 
function 

$(7,0,) = U(E,) — TS(E,) . (78) 
that depends on temperature in the usual way. We 
thus set up a one-parameter family of curves 
arranged as in Figure 10, the members of which 
respectively at successive temperatures make up the 
equilibrium free energy curve of the system. From 
this point of view the terminal structures zero and 
one have essentially the same status as £,, only the 
latter appears as the equilibrium phase at just one 
specific temperature 7}. 


To calculate 7,, we argue as follows. The phase 
¢, formally undergoes an infinitesimal first order 
transition to ¢, + df, in the neighbourhood of 7,, 
whence 


W(T;,5; + d¢,) H(7,,¢1), 
: dU 
1.@., dé, dé, = 1 aC, = @. are (19) 
Since U(,) and S(¢,) are known, Equation 19 yields 
10 
o 
208 
° 
06 
¢& 
= 
0-4 
oi] 
© 
S02 











Temperature 


Figure 9. Temperature dependence of long 
range order parameter 


T, as a function of £,, and hence also the equili- 
brium free energy of the system ¢(7;,5,) at 7;. In 
the usual formulation of the theory, emphasis is 
placed on a particular temperature 9, and the 
minimum free energy at © selected from the set 
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(9,1), . . . 6(9,0,), . . . (0,0) 
by the equation 
dU dS 
Fi 2) c = 1o)] _¥ 0, 
$( »6) d¢ d¢ 


. 


which provides the equilibrium ¢£ in terms of 0. 
Evidently this formulation is quite equivalent to the 
preceding, the former being perhaps philosophically 
superior in enabling the second order transition to 
be identified as a continuous series of infinitesimal 
first order transitions. It may be remarked that the 
latent heat of transformation at any stage is U’(¢)d¢, 
so that the total latent heat is 





¢=0 = T, 
“dU “dU 
a6 = i 
ja * ~jar* 
t=1 Te =T,—4T 


This is absorbed over a finite temperature interval, 
and hence identified as a catastrophic increase in 
the specific heat. 

Notwithstanding that the product PV may be 
neglected in comparison with U and TS, first order 
transformations are always accompanied by a 
macroscopicallly discernible change in volume. This 
feature, indeed, sometimes provides the most direct 
indication that a transformation has taken place. 
By contrast, second order transformations do not 
involve a volume change since the atoms merely 
undergo redistribution over a common lattice. The 
use of the term second order is seen to be entirely 
consistent with the greater subtlety inherent in the 
latter situation. 


The Activation Energy Barrier 
Most irreversible processes inevitably bring into 
play forces which, to a greater or lesser extent, 
oppose the reaction. At low temperatures, and 
sometimes even at room temperatures, these may 
completely inhibit the approach to equilibrium, as 
with glass or quenched-in lattice defects. At higher 
temperatures, however, the thermal energy of the 
atoms suffices for them to surmount the ‘ activation 
energy barrier’ and thus enable the reaction to 
proceed. A useful picture of the situation is pro- 
vided by the single particle model. Instead of the 
side of the well being perfectly smooth, as envisaged 
in Figure 11, we may picture it more realistically as 
broken by local hollows and humps which impede 
descent. 

Suppose the particle rests in local equilibrium at 
a potential level ¢, and that this is separated from 
the next level of local equilibrium at d6+dd¢ by a 
potential hump of height A. Although potentia! 
energy can be reduced by the transition from ¢ to 
¢+dd¢, no direct mechanism exists for overcoming 
the barrier A separating the two configurations 
However, an indirect mechanism is provided by 
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thermal vibrations. In effect, these allow the particle 
te occupy a series of harmonic oscillator levels 0, 
€, €,..., Witha relative chance e~4/A7 of acquiring 
sufficient thermal energy to jump the barrier. 
Sooner or later this latter contingency will be 
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Figure 10. (a) Free energy curves (1), 

$(¢,), $0), referring to structure of perfect 

order, partial order €, and complete disorder, 

respectively; (b) Close-up view of situation 

in neighbourhood of temperature T,, where 

$(¢,) crosses $(€, —d€,) and is in turn crossed 

by $(f + dg) 

realized, and the particle descends to the next 
configuration ¢+d¢, so continuing until final 
equilibrium has been achieved. 


It may be noted that e-*/‘T approaches unity 
when the temperature T gets large, in other words 
each level « tends to have an equal chance of being 
occupied, with the result that the higher energy 
levels are favoured at the expense of the lower. 
Consequently, the particle acquires a progressively 
increasing chance of surmounting the barrier. This 
conclusion forms the quantitative basis for under- 
standing the dominant influence of temperature on 
the rate of reaction. 


The preceding ideas have an immediate applica- 
tion to the theory of diffusion, particularly the 
diffusion of solute atoms through a crystal. Any 
solute atom lies within a potential well of height g, 
cefined by the periodic field of the crystal (Figure 


12), with a relative chance e—4¢*T of jumping into 
a neighbouring position. On multiplying this chance 
by », the frequency of atomic vibrations, we arrive 
at 4ve—4/AT for the number of successful jumps per 
second in a particular direction, bearing in mind 
that there are three equally likely independent 
directions of vibration. Typically, g is of the order 
one electron-volt, KT of order one fortieth of an 
electron-volt at room temperature, and »v of order 
10'° sec-!. If now there are n solute atoms per unit 
volume, concentrated in lattice planes of spacing 
‘a’ in the diffusion direction, this implies na solute 
atoms per unit area of lattice plane: 4na of these 
will, on the average, attempt to jump forward, of 
—_ tna. Jve—VkT ..+.(20) 
per unit time will be successful. The diffusion co- 
efficient is usually defined to be 

D = Jva® e~ 47, ies al) 
in terms of which (20) may be written nD/a. By 
measuring the variation of the diffusion coefficient 
D with temperature, or rather log D against 1/7, 
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(b) 
Figure 11. (a) Configuration of local 
equilibrium ¢ separated from neighbour- 
ing configuration ¢ + dd by a local barrier 
A; (b) Close-up view of situation at ¢. 

By virtue of the thermal levels 0, €,, € 
...A... the particle has a certain 
chance of surmounting the barrier 
we may form an estimate of g and v. It will be noted 
that the diffusion coefficient depends on the partic- 

ular direction under consideration. 


We now extend the argument to allow for the 
effect of an external potential field 4 interacting 
with the solutes, e.g., as in A. H. CoTTRELL’s theory 
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of atmosphere formation around dislocations. 
Supposing ¢ to decrease slowly over a few inter- 
atomic spacings in the diffusion direction, its effect 
will be to damp the amplitude of the periodic 
crystal field in that direction, and so impose a 
statistical directionality on the otherwise random 
thermal agitation of the atoms. Quantitatively 
formulated, if the activation energy barrier dimin- 
ishes from q+¢4, to q+d¢. over an interatomic 
spacing as shown (Figure 12), the excess relative 
chance of jumping down the gradient is 


e-@ +by)/kT __ e ~(q +) /kT 
which equals 
Gi— $s oak _ a Wo -aikr «s€22) 
kT kT dx 


on bearing in mind that ¢,—¢, is equal or nearly 
equal to —a.dédx and assuming ¢ to be small in 
comparison with kT. 


Multiplying (22) by the same factors as previously, 
we arrive at 

nD dd 

kT dx 
for the net number of solutes per unit area drifting 
down the potential gradient per unit time. The 
‘ drift velocity ’ is readily seen to be given by 
_D dd 

kT’ dx 
a formula originally deduced by Einstein in his 
theory of Brownian motion, though with a mobility 
factor arising from molecular collisions instead of 
from a periodic lattice potential. It should be noted 
that D/kT approaches zero both when T is either 
large or small: at high temperatures the solutes 
readily jump the barrier in either direction, regard- 
less of the perturbing potential ¢, so that v has zero 
value on the average; at low temperatures the 
solute atom remains trapped in its potential well, 
so that again v is zero. At intermediate tempera- 
tures the combination of the two opposing factors, 
namely the energy barrier and thermal fluctuations, 
produces the effect of a viscous drag impeding the 
motion of the solute and imparting to it an average 
limiting velocity given by Equation 23. 

Similar arguments show that — D.dc/dx gives 
the number of solutes »er unit area drifting per 
unit time down a corcentration gradient. The 
corresponding drift velocity is 

— 2) ‘de 
ce dx’ 


v 


» (29) 


Metastability 


Cases often arise where the activation energy barrier 
is so high as to be completely immune to thermal 
fluctuations. The system then carries on existing in 
the relevant state of local equilibrium, with no 
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prospect of ever achieving the greater stabiliiy 
potentially available. This phenomenon is termed 
metastability (metastable equilibrium). The situ- 
ation in terms of the one particle model appears as 
indicated in Figure 13. Perhaps the most outstand- 
ing example of a metastable system is diamond, 
which co-exists at all temperatures with the more 
stable modification graphite (Figure 14). 

The thermodynamic status of diamond differs in 
several important respects from that of a non- 
equilibrium system such as glass. Some of the 








(a) (b) 

Figure 12. (a) Periodic crystal field; (b) 
effect of slowly varying perturbing potential 
differences can be appreciated directly by comparing 
the respective one particle models: diamond has a 
clearly defined internal structure, whereas glass has 
not; diamond remains locally stable against any 
thermal fluctuations that may occur, whereas glass 
is locally unstable against thermal fluctuations that 
lead to increasing order. Diamond, in contrast to 
glass, obeys an equation of state and has a vanishing 
entropy at O°K. If graphite did not exist, diamond 
would be regarded as the stable equilibrium state of 
crystalline carbon over the whole temperature range 
of its existence. This example makes it clear that 
stable and metastable equilibrium are purely rela- 
tive terms, their use depending on the energy level 
of the reactions under consideration: no genuine 
distinction can be drawn between the two, and all 
the principles of thermodynamics, e.g., the third 
law, have the same validity for the one as for the 
other. 


Metastable systems abound in nature, indeed they 
comprise most of the systems of actual experience. 
Diamond has already been mentioned. Martensite, 
the main hardening constituent of steels, is meta- 
stable relative to a mixture of ferrite and cementite. 
Most metals are metastable, i.e., the mixture metal 
plus oxygen is metastable relative to the metallic 
oxide. Supercooled liquids are metastable relative 
to the crystalline phase, and supersaturated vapours 
are metastable relative to the liquid or crystal. A 
mixture of hydrogen and oxygen is metastable 
relative to water. On the microscopic scale, metals 
owe their mechanical properties to metastable lattice 
defects known as dislocations. Organic molecules 


are weakly metastable, thus enabling energy to be 
trapped and made available for use by the organism 
when required. As regards the kinetics of chemical 
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THE THERMODYNAMIC 


change, Sir C. N. HINSHELWoop!® has remarked that 
‘the real unities underlying the interpretation of 
chemical changes are seldom open and apparent but 
tangled in a confusion which here and there might 
well appear inextricable, but did not a broader view 
reveal that through the tangle run certain threads 
bright enough to show a plan. The strongest and 
brightest of these threads is the idea of the activa- 
tion energy’. 


Whence arises the activation energy barrier which 
maintains the various metastable systems enumer- 
ated? An instructive answer is provided by con- 
sidering the mechanism of phase nucleation in 
condensed systems, e.g., that of a crystal from 
supercooled liquid. According to the generally 
accepted picture, we envisage embryos of crystal 
phase created by random statistical fluctuations of 
atomic structure within the liquid. These embryos 
are subject to two energy factors, the first is the 
driving energy for the transformation from liquid to 
crystal, i.e., the free energy difference, 4, per unit 
volume between crystal and liquid which is markedly 
temperature dependent: at the equilibrium temp- 
erature 7, between liquid and crystal, ¢ has the 
value zero; above 7, it is positive, and below 7, 
it is negative. The second is the surface energy o 
per unit area created by the disregistery of structure 
at the boundary of the embryo: this is always 
positive and largely independent of 7. Accordingly, 
the energy required to create an embryo of radius 
r at temperature T is given by 

a=a(r,T) = or + or* 
apart from unimportant numerical factors. At any 
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Figure 13. Metastable equilibrium. The 

particle rests in local equilibrium at $,, with 

no chance of surmounting the barrier A 

separating it from the configuration ¢, of 
greater stability 


temperature T above 7,, this function increases 
steadily with r so that the embryo is always unstable. 
For any T below 7,, the function at first increases 
but passes through a maximum 


A= A(T) a*/¢? . (24) 
and at a critical radius 
R R(T) o/¢, . -€25) 


as sketched in Figure 15. This maximum defines the 
activation energy which must be supplied to the 
embryo before it can qualify as a stable nucleus of 
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crystal phase. It is the presence of the barrier A 
which ensures that supercooled liquid is everywhere 
in local internal equilibrium. 


Bearing in mind that || increases with the degree 
of supercooling, we can see from (24) and (25) 





Graphite 








T 
Figure 14. Free energy curves of graphite 
and diamond (schematic) 

that A and R both decrease as the temperature de- 
creases. The equilibrium number of nuclei at 
temperature T is provided by the Boltzmann factor 
e—A/kT; arguments can be advanced to show that 
A/T decreases, and hence the Boltzmann factor 
increases as 7 diminishes. Accordingly, as the 
temperature is lowered below 7., two possibilities 
will arise. At least one stable nucleus is created 
almost at once by virtue of a fluctuation and 
having reached the critical size then grows 
irreversibly until the phase change has been com- 
pleted. Alternatively, no stable nucleus is created 
until a temperature 7’ well below 7,: within this 
temperature interval, the liquid is said to be super- 
cooled, having then the status of a metastable 
system relative to the crystalline phase. 


The above analysis envisages nucleation as a 
random process, having an equal chance of taking 
place anywhere throughout the medium. In prac- 
tice, however, it is usually triggered off at impurities, 
at defects or at boundaries: these provide concentra- 
tions of available energy, which readily enable the 
requisite local fluctuations in structure to be effected. 
Indeed, without such preferred sites, supercooling 
and supersaturation would become much more 
pronounced phenomena than is actually the case. 
For instance water usually freezes at 0°C, but by 
careful control of size and impurities’ it can be 
supercooled through as much as 40°C. According 
to Frank’s theory of crystal growth from the vapour, 
crystallization only takes place at the observed rate 
because suitable nucleation centres are provided by 
dislocations. Gas reactions are usually initiated at 
cracks and cavities in the walls of the container. 
The two mechanisms of nucleation are termed 
homogeneous and heterogeneous respectively. 

For reactions within solids, the behaviour of a 
nucleus is controlled by a third factor, the elastic 
strain energy, besides the other two. This arises 


because the product phase usually has a different 
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specific volume from that of the matrix, thus 
generating accommodation stresses in the latter. 
Strain energy, unlike the surface energy factor, is 
proportional to the volume, and could therefore 
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Figure 15. T,>Te>T,>T3 








remain important even when the nucleus has grown 
into a precipitate of macroscopic dimensions. This 
applies particularly to martensite plates, owing to the 
coherency requirements at the boundary between 
matrix and product. The thermodynamics of 
martensite systems provides a fascinating example of 
interplay between chemical and mechanical energy, 
resulting in such phenomena as rubber-like elasticity 
and thermo-elastic equilibrium. More is likely to 
be heard about these aspects in the future, but they 
cannot be pursued further within the limits of the 
present article. 


The Thermodynamic Field 


We now introduce the concept of the thermo- 
dynamic field. This may be pictured as a medium 
where the thermodynamic variables, such as pressure 
and temperature, depend on position in a prescribed 
way. Accordingly, 
a P(x, y, z) and T = T(x, y, z) 
where x, y, z denote the coordinates of a point. A 
small system of given phase structure thus has a 
free energy 
d = o(P,T) = d[P(x, y, z), T(x, y, z)] 

P(x, y, 2), . (26) 
which varies with position, and will therefore tend 
to drift through the medium towards the nearest 
local minimum of ¢. By analogy with the drift of 
solute atoms through a crystal, we write 

D 
i iT Ve. 


For most macroscopic systems, either D or y¢ is 
so small that the system remains fixed at some 
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point, in stable or metastable equilibrium relative to 
the variables P and T operative at the point, but 
formally in non-equilibrium relative to the thermo- 
dynamic field as a whole. 


An interesting example of macroscopic finite driit 
occurs with servo-mechanism devices. Here the 
system contains a store of potential energy coupled 
to the local field conditions, and released according 
to some specified law, for example, so as to propel 
the system in the direction of decreasing tempera- 
ture. Identifying ¢ of (26) as the stored potential 
energy, the servo-mechanism may be regarded as a 
system moving in a resisting medium under the 
influence of a field of force d. Much of the theory 
of servo-mechanisms (cybernetics) is taken up 
with the way information is signalled from environ- 
ment to system, and with how the system responds. 
An ordinary dynamical system differs from a servo- 
mechanism only in the fact that it is coupled directly 
to the field, e.g., gravity, in a formal way that 
cannot be understood on the basis of any physical 
picture. Elementary biological organisms may be 
identified as servo-mechanisms in the first approxi- 
mation, the energy store being chemical rather than 
electro-mechanical. Their behaviour has _ been 
variously described as goal-seeking, purposeful 
teleological, homeostatic. However, within the 
framework of thermodynamic principles, nothing 
more is involved than drift down a free energy 
gradient. 
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SURVEY 


Recording Rail Track Behaviour 
Under Dynamic Loading 


IT is axiomatic that before passenger trains can be 
driven at high speed the track on which they are 
to run must be laid accurately and maintained in 
that condition. Early track recording gear which 
was of limited value has now been superseded, so 
far as British Railways are concerned, by a new 
coach recently completed by Elliott Brothers 
(London) Ltd. The coach has been under develop- 
ment for about two years, though a design study was 
commenced earlier to facilitate preparation of a 
specification. 


The new coach provides means for assessing 
track conditions under a loaded vehicle while it is 
moving at speed along the track; this enables track 
to be surveyed while normal timetables are in 
operation, and British Railways stipulated that 
measurements should be made in such a way as to 
make them independent of the movement of the 





Figure 1. The Elliott track recording coach 
—the photograph illustrates the sensing shoe 
along the running rail 


vehicle. The coach has been designed as a four- 
wheeled self-propelled vehicle, operable from either 
end, to record at any speed up to 30 miles per hour 
and when not recording to travel at speeds up to 
55 miles per hour. One axle only is driven; the 
other is used to provide mechanical drive to the 
distance-measuring apparatus. 


Continuous records (traces on photosensitive 
paper) of superelevation (cant), curvature and 
gauge are produced, and additional equipment is 
provided so that a record of an earlier journey by 
the coach over the same track can be run through in 
synchronism. This enables engineers to see immedi- 
ately whether track adjustments which may have 
been made in the interim have had the desired 
effect. Compensating adjustments can be made on 
the recording apparatus to take account of the 
continuous rotation of the earth (every half-hour) 
or the motion of the coach around curves. This is 
necessary because the datum for cant measurement 
is provided by a high speed gyroscope indepen- 
dently of curvature or gradient: this is compared 
with the tilting of one axle. Curvature and gauge 
are obtained by measuring the movements of a 
system of probes which contact the inside edge of 
the head of the rail. 


These probes (Figure /) are designed to negotiate 
track fittings at the maximum recording speed. The 
rubbing faces of the sensing shoes are of Stellite, 
and it is claimed that wear is so small that no 
allowance need be made for this even on tests 
exceeding 1000 miles in length. It should be noted 
that the cone to which the running wheels of this 
coach are turned is 1:100. Normal wheels would 
introduce too much variability into the record of 
cant. Normal loading is about 12 ton/axle, and 
the coach carries its own generating equipment to 
supply the gyroscopes and synchros at 400 c/s, 
three phase. Living accommodation for the crew 
is also provided. 


Industrial Application of Fluid Flow 


The British Hydromechanics Research Association 
recently held Open Days at their establishment at 
Harlow, Essex. This Association is non-profit 
making, being controlled and financed by members 
with the help of a proportionate government grant. 
Any industrial applications of fluid flow within fixed 
boundaries is of interest to the Association and 
problems investigated include: mine ventilation; 
centrifugal pump design; water hammer in pipe- 
lines; seals, packings and joints; cavitation; 
hydraulic transport of solids; friction losses in pipe 
fittings; and the hydraulic design of civil engineer- 
ing structures. In fact, during work on the last- 
named problem, it has proved possible so to modify 
channel design and control valve equipment as to 
make large savings of capital on initial outlay and 
to reduce running costs substantially by eliminating 
turbulence, removing potential vortices and reducing 
frictional and head losses in various ways. 
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BOOK REVIEWS 


A History of Western Technology 
F. KLEMM 
(401 pp; 8} in. by 5% in.) 
London: George Allen & Unwin. 32s 


THERE are various ways_in which one might attempt to 
write a history of technology. For instance, one could 
describe the principles involved in the development of 
machines and techniques; one might discuss the in- 
dividual inventors, and innovators of techniques, giving 
their biographies and an assessment of the importance 
of their work; or one could, as Mumford has done, 
write a history of society and its dependence on the 
development of machines. 


In writing A History of Western Technology, Dr 
KLEMM has done none of these things. In fact, his book 
may perhaps be more properly described as an anthology 
of carefully selected contemporary writings connected by 
commentary in the text, rather than as a history. Dr 
Klemm is a man of very wide erudition, and his extracts 
range from those taken from technical treatises such as 
De Re Metallica (from which he gives a good deal of 
the most interesting introduction) to quotations from 
contemporary political and religious writings. His in- 
tention appears to be not only to describe the progress 
of technology but also to outline the background condi- 
tions under which this progress occurred. In this aim 


he succeeds admirably; and at the same time he gives 


the reader access to material and documents which are 
not always readily available. It would not, one feels, 
be going too far to say that this method is ideal if one 
wishes to outline the progress of technology (or science) 
and to deal with the subject in a humanistic way. 


There is, however, one slight criticism of this other- 
wise admirable book that one might venture. It is a 
relatively short book, but in it the author attempts to 
cover the whole of man’s technological development 
from Grecian times to the detonation of the first atomic 
bomb. Somewhat naturally, therefore, the quality of 
the treatment varies. The earlier periods up to and 
including the Renaissance are very good; but the later 
periods tend to be treated rather sketchily. One feels 
that Dr Klemm has tried to cover too much ground 
and that it would have been better to have ended the 
book at about 1900 (at the end of what Mumford calls 
the ‘ Palaeotechnic Phase ’) rather than to have ventured 
into the present phase of extremely rapid technological 
development. 


This is a history of Western, and not British, techno- 
logy. Inevitably, therefore, there are omissions and 
inclusions which, from our insular viewpoint, we may 
find surprising. Some of these are obviously salutory as, 
for instance, the full discussion of Papin’s contribution 
to the development of the steam engine which is often 
overshadowed here by the work of Savery and New- 
comen. Less excusable in this context, however, is the 


absence of any reference to Robert Hooke who was no 
only Papin’s guide and mentor for many years but who 
also made many very important contributions to 
technology. 

Nevertheless, this book can be thoroughly recom- 
mended. It is well written, and the translation is 
extremely good. The plates and diagrams are well 
chosen and extremely well produced. The list of sources, 
the bibliography and the index are excellent; and as 
was said above, the texts quoted have been very carefully 
chosen and include translations from many which are 
difficult of access in this country. The book fills a gap 
which many of us are conscious of on our shelves; it 
provides information and diagrams which few of us have 
readily available; but above all it is a notable addition 
to the all too sparse literature on this most important 
subject. H. D. TuRNER 


Meson Physics 
R. W. MARSHAK 
(viii+-378 pp; 8 in. by 5 in.) 
New York: Dover Publications; 
London: Constable. 16s (paper covered) 


Tuis book is a republication in cheaper form of a volume 
first published in 1952 by the McGraw Hill Book Com- 
pany. In 1952 Professor Marshak’s book was welcomed 
as one of the first critical surveys of the rapidly advancing 
subject of meson physics. It must now be judged against 
more recent books, for example, Professor Jackson’s 
The Physics of Elementary Particles published in 1958. 
Comparison makes one doubt the wisdom of republica- 
tion in an unabridged and unaltered form since progress 
in this field has been so great. When Meson Physics was 
originally written the contribution of cosmic rays to pion 
and muon physics was considerable, and to the new 
particles, K-mesons and hyperons, all important and 
decisive. But with the coming into operation of the great 
proton synchrotrons at Brookhaven and Berkeley in 
1953 and 1954, and the development of new techniques 
such as the diffusion chamber and the bubble chamber, 
the position has radically changed. New particles have 
been discovered and a wealth of factual information has 
been obtained. 

Advances on the theoretical side since 1952 have been 
no less impressive. Of especial significance has been the 
discovery of the non-conservation of parity in weak 
reactions with its important consequences for the theory 
of the neutrino and meson decay processes. Mention 
should also be made of associated production of the 
* strange’ particles, and the Gell—Mann-Nishijima 
scheme. 

Re-reading Professor Marshak’s book after a lapse of 
six years one is impressed with his success in capturing 
the drama of the challenging enterprise which is modern 
particle physics. The drama is still being enacted but 
those who know of recent progress will wish for Acts I 
and II and not just Act I. G. D. ROCHESTER 


See also page 273 
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